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SUMMARY j 

Recent  interest  in  a  wide  variety  of  low  Reynolds  numher  configurations  has  focused  attention  on  the 
design  and  evaluation  of  efficient  airfoil  sections  at  chord  Reynolds  numbers  from  about  100,000  to  about 
1,000,000.  These  configurations  include  remotely  piloted  vehicles  (RPV'sl  at  high  altitudes,  sailplanes, 
ultra-light  man-carrying/man-powered  aircraft,  mini-RPV  s  at  low  altitudes  and  wind  turhines/propellers.  A 
study  is  presented  of  the  present  status  and  future  possibility  of  airfoil  design  and  evaluation  at 
subcritical  speeds  to  meet  the  needs  for  these  applications. 


Although  the  design  and  evaluation  techniques  for  airfoil  sections  above  chord  Reynolds  numhers  of 
500,000  is  reasonably  well  developed,  serious  problems  related  to  boundary  layer  separations  and  transition 
have  been  encountered  below  R^  =  500,000.  Presently  available  design  and  analysis  methods  need  to  improve 
their  criteria  for  laminar  separation,  transition,  and  turbulent  separation.  Improved  mathematical  models 
of  these  complex  phenomena  require  additional, very  careful  experimental  studies.  Rpcause  of  the  sensitivity 
of  the  low  Reynolds  number  airfoil  boundary  layer  to  free  stream  and  surface-  generated  dusturhances , 
definitive  experiments  are  very  difficult.  Also  the  physical  quantities  measured  (i.e.,  pressure  difference 
and  drag  forces  etc.)  are  very  small  and  the  accuracy  of  such  measurements  depenrfcon  the  method  used.  The 
results  from  numerous  experimental  studies  are  presented  to  illustrate  the  type  of  difficulties  encountered. 

1.  INTRODUCTION 

Recently  attention  has  turned  toward  low  Reynolds  number  aerodynamics  at  subcritical  speeds  in  an 
effort  to  obtain  better  performance  for  both  military  and  civilian  applications.  These  applications  includp 
remotely  piloted  vehicles  (RPVs)  at  high  altitudes,  sailplanes,  ultra-light  man-carrying/man-powered 
aircraft,  mini-RPVs  at  low  altitude,  and  wind  turbines. 


Since  the  airfoil  section  forms  the  basic  element  in  the  design  of  a  wing,  propellpr  or  conventional 
wind  turbine,  it  has  been  the  focus  of  most  of  the  attention  in  this  area1"3,  the  performance  of  an  airfoil 
section  is  critically  dependent  upon  the  character  of  the  viscous  boundary  laypr  (i.e.,  laminar, 
transitional  or  turbulent).  The  Reynolds  number,  which  is  the  ratio  of  inertial  to  viscous  forces  or 
characteristic  length  times  velocity  divided  by  kinematic  viscosity,  is  usually  used  to  scale  vehicles.  Thp 
characteristic  length  in  this  case  is  the  airfoil  chord.  Other  things  equal,  the  character  of  the  boundary 
layer  has  been  found  to  he  dependent  on  the  magnitude  of  the  Reynolds  number.  A  broad  perspective  on  the 
range  of  chord  Reynolds  numbers  versus  flight  velocity  and  Mach  number  for  a  variety  of  natural  and  man-madp 
flying  objects  can  be  obtained  from  Figure  1. 

Although  the  designer  of  a  large  transport  aircraft  might  consider  a  chord  Reynolds  numhpr  of  lh6  to  hp 
low,  the  designer  of  a  high  altitude  RPV  or  wind  turbine  would  take  the  opposite  point  of  view.  Almost  all 
natural  flying-objects  fall  in  the  Reynolds  number  range  below  in6.  Furthermore  it  should  be  remembered 
that  most  large  high  speed  vehicles  or  at  least  some  components  of  these  vehicles  (i.e.,  control  surfaces 
and  lift  or  drag  augmentation  devices)  operate  at  much  lower  Reynolds  numhprs  durinq  takp-off  and  landinq. 
It  should  also  be  mentioned  that  much  of  the  wind  tunnel  testing  by  large  aircraft  manufacturers  is  done  at 
relatively  low  Reynolds  numbers.  At  high  altitudes  aircraft  gas  turbine  engine  fan,  compressor,  and  turhinp 
blades  with  their  small  chords  encounter  Reynolds  numbers  considerably  below  106.  Fypn  the  Space  Shuttle 
encounters  Reynolds  numbers  as  low  as  10^  at  M  =  21  during  reentry.  It  is  clear  that  some  definition  of 
what  is  meant  by  low  Reynolds  number  aerodynamics  is  needed.  In  the  present  context,  low  Reynolds  numhers 
are  considered  to  be  those  below  about  in6  since  the  effects  of  laminar  separation  and  transition  to 
turbulence  in  the  airfoil  boundary  layer  have  a  great  and  sometimes  unpredictable  influence  on  overall 
vehicle  performance.  The  applications  included  under  the  titlp  low  Reynolds  Number  Vehicles  require 
efficient  airfoil  sections  in  the  chord  Reynolds  number  range  from  about  in6  to  ID6  as  indicated  in  Figure 
1. 


1.1  Applications  of  Interest 


A  large  number  of  applications  have  been  proposed  for  high  altitude  RPVs,  often  rpferrrd  to  as  hiqh 
altitude  aircraft  platforms  (HAAPs)6"1  .  Braves11  summarized  thesp  proposed  haap  applications  as  follows: 
1)  Military  -Communications  Relay,  Ballistic  Missile  Farly  Warning,  Aircraft  Trackinq,  Weather  Monitorinq, 
Ocean  Surveillance,  Rattlefield  Tactical  Intelligence,  and  Nuclear  Explosion  Cloud  Samplinq,  2)  Scientific 
-Astronomical  Observations,  Atmospheric  Research,  and  Oceanographic  Research,  1)  Civil  -200  Mile  Fishery 
Enforcement,  Border  Patrol  Surveillance,  Water  Pollution  Monitoring,  Resource  Management,  MHF  TV  Broadcasts, 
National  TV  Ulstrlbutlon,  Ice  Surveying/Mapping  of  Waterways,  and  Emergency  Response  Communications.  Since 
most  of  these  applications  would  require  that  the  vehicle  fly  continuously  without  refueling,  the 
feasibility  of  using  solar-voltaic,  microwave  and  nuclear  propulsion  systems  have  been  explored  6_1'. 


The  modern  sailplane  represents  man's  greatest  triumph  in  aerodynamic  performance  and  efficiency. 
Although  soaring  began  In  the  late  1800's,  the  vast  majority  of  improvements  in  this  technology  have  been 
made  recently  in  Europe,  especially  in  Germany  Sailplane  nevelopments  have  provided  the  stimulus  for 
the  practical  application  of  laminar  flow  aerodynamics. 


"Department  of  Aerospace  and  Mechanical  Engineering 


The  recent  ultra-light  man-carrying/man-powered  aircraft,  developed  largely  in  the  United  States,  owe 
their  success  mostly  to  lightweight  materials  and,  to  a  lesser  degree,  to  aerodynamics.  The  addition  of  a 
small  gasoline  engine  to  a  hang  glider  to  produce  an  ultra-light  aircraft  was  accomplished  as  recently  as 
the  early  1970‘s  23.  The  first  entry  was  the  so-called  "hack  pack"  engine.  A  in  hp.  Go-Kart  engine  driving 
a  caged  propeller  was  simply  strapped  to  the  back  of  the  flyer  of  the  Rogallo  "kite."  nue  to  its  weight, 
discomfort  and  poor  performance,  it  was  never  really  accepted.  Soon  after  the  "back  pack"  Rogallo  concept, 
super-lightweight  versions  of  more  conventional  airplane-like  hang  gliders  with  control  systems  were 
developed.  In  1975  a  go-kart  engine  was  mounted  to  a  biplane  hang  glider  of  this  type  by  Moody  20  and  the 
ultra-light  airplane  was  born.  In  pursuit  of  the  “Kremer  Prize",  MacCready  and  his  associates  designed  and 
successfully  flew  one  of  the  earliest  man-powered  (i.e.,  muscle-driven)  airplane.  On  August  ?3,  1977,  the 
Gossamer  Condor  flew  the  figure-eight  Kremer  course  in  7  1/2  minutes.  The  bicycle-type  propulsion  mechanism 
driven  by  B.  Allen  produced  about  0.4  hp  at  cruise.  This  historic  man-powered  flight  was  followed  by  the 
two  hour  and  49  minute  flight  of  the  Gossamer  Albatross  across  the  English  Channel  in  1979  and  later  by  the 
flight  of  the  solar-powered  Solar  Challenger  21-23. 

Mini-RPVs  flying  at  low  altitudes  may  be  used  as  target  vehicles,  decoys  and  for  battlefield 
surveillance  24-28.  These  vehicles  are  similar  in  size  to  highly  efficient  model  airplanes  usually  wi*-h 
electronic  payloads  weighing  less  than  100  kg. 

The  growing  interest  in  wind  energy  conversion  using  wind  turbines  provides  still  another  important 
application.  Although  a  wide  variety  of  windwheel  devices  have  been  used  foi  many  years  20,  recent,  interest, 
has  been  concentrated  on  the  vertical  axis  Oarrieus  rotor  and  on  the  horizontal  propeller  type  rotor. 

1.2  Design  Aims 

The  design  requirements  cover  a  wide  range  when  one  considers  the  diversity  of  applications  from  RPVs 
to  wind  turbines.  Almost  all  low  Reynolds  number  vehicles  share  the  ultimate  goal  of  maximum  aerodynamic 
efficiency.  An  aerodynamic  efficiency  map  originally  drawn  by  McMasters^q  is  shown  in  Fiqure  2.  Areas 
assumed  to  cover  the  ranges  of  operation  of  high  altitude  aircraft  platforms  and  low  altitude  mini-RPVs  have 
been  added  to  this  figure.  Since  data  of  this  type  is  generally  unavailable  for  wind  turbines,  this  item 
has  not  been  included.  A  useful  comparison  of  vehicle  weight  for  a  qiven  winq/rotor  area  obtained  from 
McMasters  30  is  presented  in  Figure  3.  In  the  region  of  interest  in  this  study,  vehicle  mass  greater  than 
about  50  kg,  the  Ma  S^/2  line  captures  the  major  trend.  Two  of  the  larqpst  existing  wind  turhines,  a 
recently  proposed  solar  powered  HAAP  and  an  estimated  region  for  mini-RPVs  are  also  included  on  Figure  3. 

Requirements  for  a  typical  low  altitude  mini-RPV,  for  example,  include  long  flight  duration  (i.e.,  hiqh 
value  of  C|3/2/Cq)  at  speeds  between  10  and  25  meters/sec  (i.e.,  chord  Reynolds  numbers  from  about  100,000 
to  500,000)  and  altitudes  from  30  to  300  meters,  light  weight,  and  all-weather  capabi  1  i ties-33-34  (j.p,, 
precipitation,  wind  shears  and  unsteady  wind).  Minimum  wing  area  for  ease  of  packaging  and  pre-launch 
handling  is  also  important.  High  altitude  RPVs  have  somewhat  higher  speed,  chord  Reynolds  number  and 
endurance  requirements;  however,  their  cruise  altitudes  are  usually  frpe  of  the  extremes  of  precipitation, 
wind  shear  and  unsteady  wind  found  near  the  earth's  surface.  They  must,  of  course  survive  the  initial 
journey  from  launch  near  sea  level  to  their  high  altitude  station.  Wind  turbine  hlades  also  require  high 
aerodynamic  efficiency  and  all-weather  capabilities.  The  need  for  efficient  low  Reynolds  number  airfoils 
which  are  not  overly  sensitive  to  wind  shear,  gusts,  and  the  roughness  produced  by  precipitation  is  common 
to  most  mini-RPV  vehicles.  Furthermore,  confidence  that  the  operational  vehicle  will  pprform  as  designed  is 
important  in  all  applications. 

1.3  Flow  Problems 

Although  design  methods  developed  over  the  past  20  years  produce  efficient  airfoils  for  chord  Rpynolds 
numbers  greater  than  about  500,000  these  methods  are  generally  inadequate  for  chord  Reynolds  numbers  below 
500,000.  In  relation  to  the  airfoil  boundary  layer,  important  areas  of  concern  are  the  separated  reqions 
which  occur  near  the  leading  and/or  trailing  edges  and  transition  from  laminar  to  turbulent  flow  ’h.  It  is 
well  known  that  separation  and  transition  are  highly  sensitive  to  Reynolds  number,  pressure  qradient,  and 
the  disturbance  environment  36.  Transition  and  separation  play  a  critical  role  in  determining  the 
development  of  the  boundary  layer  which,  in  turn,  affects  the  overall  performance  of  the  ai rfoi 1 ^8-33.  The 
aerodynamic  characteristics  of  the  wing  and  other  components  in  turn  affect  the  static,  dynamic  and 
aeroelastic  stability  of  the  entire  vehicle.  Therefore  the  successful  management  of  the  very  sensitive 
boundary  layer  for  a  particular  low  Reynolds  number  vehicle  design  is  critical. 

At  high  Reynolds  numbers,  laminar  floe  rarely  persists  very  far  downstream  of  the  leadinq  edge  as 
illustrated  in  Figure  4.  In  this  flow  regime,  the  wing  is  usually  free  of  the  relatively  large  laminar 
separation  regions  which  occur  at  low  Reynolds  numbers  as  illustrated  in  Figures  5-7.  Our  incomplete 
physical  understanding  of  the  transition  process,  which  controls  the  location  of  separation  and 
reattachment,  is  the  major  shortcoming  in  designing  this  type  of  airfoil.  At  low  Reynolds  numbers  the 
transition  process  is  much  more  sensitive  to  free  stream  disturbances  and  aberrations  of  the  airfoil 
geometry  37. 

Laminar  separation  bubbles  occur  on  the  upper  surface  of  most  airfoils  at  low  Reynolds  numbers.  These 
bubbles  become  larger  as  the  Reynolds  number  decreases  usually  resulting  in  a  rapid  deterioration  in 
performance,  i.e.,  substantial  decreases  in  (L/0).  In  principle  the  laminar  separation  huhhle  and 
transition  can  be  artificially  controlled  by  adding  the  proper  type  of  disturbance  at  the  proper  location  on 
the  airfoil.  Wires,  tape  strips,  grooves,  steps,  grit,  or  bleed-through  holes  in  the  airfoil  surface  have 
all  been  used  to  have  a  positive  Influence  on  the  boundary  ’ayer  in  this  critical  Reynolds  number  region. 
The  type  and  location  of  these  so-called  "turhulators”  and  their  actual  effect  on  the  airfoil  boundary  layer 
has  not  been  well  documented.  Furthermore,  the  addition  of  a  turhulator  does  not  always  Improve  the  airfoil 
performance38.  In  fact,  how  the  disturbances  produced  by  a  given  type  of  turhulator  influence  transition  is 
not  wel 1  understood. 


As  a  result  of  this  critical  boundary  layer  behavior,  several  important  questions  must  be  asked:  11 
What  is  the  free  stream  disturbance  level  and  flight  environment  for  a  given  low  Reynolds  numher 
application?  2)  If  the  flight  conditions  are  known  and  a  suitable  design  technique  were  available,  could 
the  resulting  vehicle  or  component  be  adequately  evaluated  in  a  wind  tunnel  which,  in  general,  has  a 
different  disturbance  level  and  environment  than  the  flight  condition?  3)  Is  the  hysteresis  in  aerodynamic 
forces  observed  in  low  turbulence  wind  tunnel  experiments  present  in  powered  applications  (i.e.,  do 
structural  vibrations  originating  with  the  propulsion  or  drive  system  affect  boundary  layer  transition)’  4) 
Because  the  critical  quantities  measured  in  wind  tunnel  experiments  are  very  small  what  is  the  level  of 
accuracy  necessary  to  improve  design  and  analysis  methods?. 

1.4  Scope  of  Review 

The  purpose  of  the  present  review  is  to  point  out  a  numher  of  design,  analysis  and  experimental 
problems  in  low  Reynolds  number  aerodynamics  which  are  of  major  importance  in  the  successful  desiqn  and 
operation  of  vehicles  in  this  regime.  Although  there  are  other  important  problems  related  to  the  vehicle 
structure,  power  system,  control  system,  etc.  this  study  will  concentrate  on  the  aerodynamic  problems  at 
subcritical  speeds. 

2.  BOUNDARY  LAYER  BEHAVIOR 

The  magnitude  and  direction  of  the  aerodynamic  forces  and  moments  actinq  on  a  movinq  vehicle  (e.g., 
aircraft  wing  or  wind  turbine  rotor)  are  determined  to  a  large  extent  by  the  behavior  of  the  boundary  layer. 
These  forces  are  related  to  how  rapidly  the  viscous  boundary  layer  grows  and  whether  or  not  it  separates 
from  the  vehicle's  surface.  The  .ate  of  growth  of  the  boundary  layer  and  whether  separation  occurs  are 
strongly  influenced  by  the  character  of  this  viscous  layer  (i.e.,  laminar,  transitional  and/or  turbulent). 
This  character  of  the  boundary  layer  is  a  function  of  the  vehicle  shape,  relative  roughness  of  the  surface, 
as  well  as  the  free  stream  disturbance  level  approaching  the  vehicle.  The  key  to  the  boundary  layer 
behavior  for  low  Reynolds  number  vehicles  is  whether  transition  from  laminar  to  turbulent  flow  takes  place 
in  the  attached  boundary  layer  before  laminar  separation  occurs.  Once  laminar  separation  occurs,  the 
subsequent  laminar  free  shear  layer  is  highly  unstahle  and  transition  to  turbulent  flow  takes  place  quite 
rapidly. 

Several  types  of  flow  behavior  result  at  low  Reynolds  numbers.  First  a  laminar  separation  may  occur 
near  the  leading  edge  at  high  angles  of  attack  in  which  case  the  airfoil  may  he  considered  fully  "stalled," 
or  at  small  angles  of  attack  the  laminar  boundary  layer  may  remain  attached  through  an  extendpd  favorable 
pressure  gradient  near  the  leading  edge  and  then  separate  in  an  adverse  pressure  qradient  after  the  maximum 
thickness  of  the  airfoil.  In  either  case  an  unsteady  oscillating  wake  is  formed  which  reduces  airfoil 
performance. 

A  second  type  of  boundary  layer  behavior  may  occur  which  is  probably  the  most  desirable  at  low  Reynolds 
numbers.  The  hest  airfoil  performance  is  achieved  when  the  laminar  boundary  layer  transitions  to  a 
turbulent  one  before  reaching  the  large  adverse  pressure  gradient.  The  turbulent  boundary  layer  with  the 
higher  energy  level  is  able  to  remain  attached  to  the  airfoil  through  the  adverse  pressure  qradient.  This 
"natural"  transition  (not  caused  by  the  separation  buhble)  is  accompanied  by  hiqher  lift  coefficients  and 
lower  total  drag  coefficients.  At  high  angles  of  attack,  a  trailing  edge  separation  of  the  turbulent 
boundary  layer  provides  gentle  "stalling"  characteristics. 

A  third  type  of  boundary  layer  behavior  occurs  which  may  he  considered  an  extension  of  the  laminar 
separation  case.  Instead  of  remaining  separated,  the  laminar  free  shear  layer  in  some  cases  may  reattach 
shortly  after  separation  or  the  free  shear  layer  more  often  may  becomp  turbulent,  and  the  growing  turbulent 
shear  layer  then  interacts  with  the  airfoil  surface  usually  causing  reattachmpnt .  After  reattachment,  the 
turbulent  boundary  layer  behaves  in  a  manner  similar  to  the  natural  transition  case  although  it  thickens 
more  rapidly.  Theoretical  prediction  of  the  presence  and  location  of  the  laminar  separation  buhble  and 
experimental  studies  of  che  behavior  of  separation  bubbles  have  been  the  focus  of  many  previous 
investigations.  Separation  bubbles  have  been  placed  in  two  classifications:  the  short  bubble  and  the  long 
bubble.  The  short  bubble  typically  occurs  at  high  angles  of  attack  and  <s  usually  less  than  a  few  percent 
of  the  chord  in  length  at  high  Reynolds  numbers.  However  at  low  chord  Reynolds  numbers  Tani™  observed 
short  separation  huhbles  that  were  approximately  2BT  of  the  chord  in  length.  A  short  separation  buhble  at 
high  Reynolds  numbers  has  very  little  effect  on  the  overall  theoretical  pressure  distribution  and  usually 
decreases  in  size  with  increasing  incidence.  In  comparison,  the  length  of  the  lonq  huhhle  is  of  more  than  a 
few  percent  of  chord  and  will  lengthen  as  angle  of  attack  increases  or  Reynolds  numher  decreases.  Thp 
presence  of  a  long  bubble  greatly  alters  the  pressure  distribution  from  its  theoretical  form.  A  separated 
laminar  boundary  layer  is  sometimes  thought  of  as  a  long  bubble  extpndinq  into  thp  wake  of  the  airfoil. 
Figure  8  from  References  40  and  41  illustrates  the  behavior  of  the  pressure  distribution  with  the  formation 
of  a  long  or  short  bubble. 

As  the  chord  Reynolds  number  decreases,  the  laminar  portion  of  the  free  shear  layer  in  a  short  bubble 
grows  in  length  and  the  turhulent  portion  requires  more  entrainment  to  reattach  at  a  pressure  near  the 
inviscid  pressure  value.  Eventually  the  Reynolds  number  becomes  so  low  and  the  laminar  portion  of  the  free 
shear  layer  so  long  that  the  turhulent  entrainment  process  can  no  longer  support  reattachment  near  thp 
inviscid  pressure  value.  The  velocity  peak  and  circulation  decrease  reducing  the  pressure  gradient  ovpt  thp 
bubble.  This  allows  the  turbulent  free  shear  layer  to  reattach  as  a  long  bubble  and  the  short  bubble  is 
said  to  nave  "burst"  into  a  long  bubble.  This  separation  bubble  decreases  airfoil  performance  (i.e., 
increases  pressure  drag  and  results  in  a  much  thicker  turbulent  boundary  layer  downstream  of  the  bubble)  as 
compared  to  the  natural  transition  case,  but  is  a  large  Improvement  over  the  separated  laminar  boundary 
layer  which  does  not  reattach.  The  presence  of  a  short  laminar  separation  bubble  can  he  utilized  to  improve 
airfoil  performance  at  very  low  Reynolds  numbers  hecause  it  acts  as  a  trip  and  reduces  the  possibility  of 
massive  separation  further  downstream. 


Separation  bubbles  often  have  a  dramatic  effect  on  the  stalling  characteristics  of  airfoils.  When  a 
short  bubble  is  present  on  an  airfoil  the  lift  increases  linearly  with  angle  of  attack  until  stall  occurs. 
A  large  discontinuous  drop  in  lift  usually  accompanies  the  bubble  hurst  at  high  Reynolds  numbers  as  shown  in 
Figure  9.  If  a  long  bubble  forms  on  an  airfoil,  stall  occurs  when  it  has  extended  to  the  trailing  edge. 
The  resulting  lift  curve  peak  is  fairly  flat  and  has  no  discontinuities.  At  low  Reynolds  numbers  these 
effects  plus  additional  problems  such  as  hysteresis  are  present. 

Seme  airfoils  exhibit  a  phenomenon  near  stall  in  which  the  aerodynamic  forces  developed  depend  on  the 
direction  the  angle  of  attack  was  reached.  As  the  angle  of  attack  increases,  the  lift  and  drag  forces 
increase.  At  stall,  an  abrupt  decrease  in  lift  and  increase  in  drag  occurs.  A  small  reduction  in  the  angle 
of  attack,  however,  does  not  restore  the  forces  to  their  former  values.  Instead,  the  angle  may  have  to  he 
reduced  several  degrees  before  the  lift  and  drag  suddenly  revert  to  the  values  obtained  under  conditions  of 
increasing  angle  of  attack.  This  behavior  shown  in  Figure  10  is  described  in  Reference  42  as  “high-lift" 
hysteresis  and  in  Reference  43  as  "high  Cynax"  or  "clockwise"  hysteresis  and  is  attributed  to  the 
development  and  bursting  of  a  short  bubble. 

The  reverse  situation  may  also  occur  in  which  an  abrupt  increase  in  lift  and  decrease  in  drag  takes 
place  at  high  angles  of  attack.  The  forces  do  not  revert  to  the  values  obtained  under  conditions  of 
increasing  angle  of  attack  until  a  sufficiently  low  angle  is  achieved  as  shown  in  Figure  11.  Reference  4? 
calls  this  “moderate-lift"  hysteresis  and  Reference  43  calls  it  "low  Cpax"  or  "counterclockwise" 
hysteresis.  This  type  of  hysteresis  appears  to  result  from  the  growth  of  a  long  bubble  and  its  sudden 
collapse  into  a  short  bubble. 

2.1  Historical  Background 


The  complex  mechanisms  of  separation,  transition  and  turbulence  have  been  the  focus  of  many 
investigations  over  the  past  half  century.  The  phenomenon  of  separation  bubbles  was  first  studied  by  R.M. 
Jones  in  193344.  He  noted  the  existence  of  a  separation  and  reattachment  of  the  flow  over  thin  slightly 
cambered  airfoils.  In  the  early  1940's,  F.U.  Schmitz  performed  numerous  experiments  on  model  airplanes 
operating  in  the  low  Reynolds  number  regime4'.  He  was  one  of  the  first  to  recognize  t.he  advantage  of  the 
separation  bubble  in  reducing  drag  hy  tripping  the  boundary  layer  and  thus  maintaining  an  attached  flow  over 
a  larger  portion  of  the  airfoil.  O.E.  Gault,  in  a  series  of  experiments  conducted  between  1949  and 
19554®.47,  investigated  the  regions  of  separated  laminar  flow  and  categorized  the  types  of  boundary  layer 
behavior.  P.R.  Owen  and  L.  Klanfer4®  studied  separation  bubble  bursting  on  a  thin  airfoil  and  made  an 
attempt  to  develop  a  criterion  for  bursting.  In  1956  U.  Pfenninger  studied  the  performance  of  gas  turbine 
blades  operating  in  the  Reynolds  number  range  of  30,000  to  100,000  and  concluded  that  performance  could  he 
improved  by  shortening  the  separated  laminar  region  artificially  by  means  of  disturbances  introduced  into 
the  boundary  layer4^. 


Moore  continuing  the  work  of  Owen  and  Klanfer  in  I960,  experimentally  determined  that  when  the  Reynolds 
number  based  upon  boundary  layer  displacement  thickness  was  less  than  500  at  separation,  that  a  long 
separation  bubble  was  formed.  When  this  parameter  was  greater  than  500  a  short  buhble  formed  although  it 
might  burst  as  the  airfoil  anqle  of  attack  increased®®1.  Raster  performed  an  extensive  study  of  the 
separation  bubble  in  the  1960's4*>®1 .  He  established  conditions  under  which  short  bubbles  could  burst  to 
form  long  bubbles. 


In  an  excellent  review  in  the  early  1960's,  J.W.  Ward  combined  experimental  and  theoretical  results 
into  a  methodology  for  studying  laminar  separation  and  bubble  format  ion5’'’.  It  was  generally  accepted  that  to 
improve  airfoil  performance  at  low  Reynolds  numbers  some  type  of  mechanism  was  necessary  to  cause  the 
laminar  boundary  layer  to  transition  to  a  turbulent  one.  The  formation  of  a  separation  bubble  was  one  of 
several  options.  Other  options  included  the  use  of  boundary  layer  roughness  or  trips,  increased  free  stream 
turbulence,  transpiration  (blowing  of  air  into  the  boundary  layer),  or  the  use  of  the  adverse  pressure 
gradient  itself.  The  formation  of  the  laminar  separation  bubble  and  bubble  hurstlnq  have  hepn  studied  by 
Tani39f  Horton^J^  Ntim®4,and  Roberts”. 


In  1965,  J.L.  vanlngen  conducted  a  theoretical  and  experimental  study  of  the  incompressible  boundary 
layer  and  methods  of  influencing  transition”.  More  recently,  n.  Althaus  has  studied  the  Influence  of 
roughness  Introduced  by  insect  remains  on  the  performance  of  aifoils  operating  at  low  Reynolds  numbers'*7. 

Up  until  the  past  few  years,  low  Reynolds  number  studies  of  boundary  layer  performance  (with  the 
exception  of  Schmitz  and  Pfenninger)  have  concentrated  on  performance  at  Reynolds  numbers  above  5(1(1,000. 
Unfortunately,  many  low  Reynolds  number  applications  fall  below  this  value.  O.F.  Volkers  has  studipd  thp 
lift  and  drag  performance  of  several  airfoils  between  Rc  =  60,000  ann  5(1(1, 00n®fl.  Arena  and  Mueller®q  and 
Mueller  and  Rurns®^  made  extensive  use  of  flow  visualization  to  study  ooundary  layer  performance  below  Rc  = 
500,000.  Most  recently  Conigl la-o®*  and  Mueller, et  al.®7  conducted  studies  of  thp  airfoil  section  used  on 
the  first  man-powered  aircraft.  Jansen®^  completed  a  thorough  study  of  the  separation  buhble  on  the  NAPA 
663-OI 8  airfoil.  A  detailed  study  of  the  boundary  layer  characteristics  of  the  Miley  airfoil  between  Rc  = 
70,000  and  600,000  has  been  completed  by  Pohlen4®.  Numerical  studies  of  low  Reynolds  number  airfoil  flow 
problems  have  recently  been  performed  for  both  steady  and  unsteady  flows®®"®  .  AH  of  these  previous 
studies  have  formed  a  basis  for  the  study  of  airfoil  performance  at  low  Reynolds  numbers  which  should  help 
in  the  design  of  airfoils  for  this  flight  regime. 

2.2  Laminar  Flow 


Laminar  flows,  whether  steady  or  unsteady,  two  or  three-dimensional,  may  be  determined  by  a  variety  of 
analytical  and  numerical  techniques.  Entire  laminar  flows  including  regions  of  separation  have  been 
accurately  determined  using  several  different  numerical  methods  of  solution  of  the  discretized  Navier-Rtokes 
equations  with  the  appropriate  boundary  conditions.  Using  modern  grid  generation  techniques7®1  for  arbitrary 
shapes,  approximate  solutions  of  this  type  are  becoming  commonplace.  Attached  laminar  boundary  layers  may 
be  calculated  with  a  much  smaller  expenditure  of  programming  and  computer  time  hy  simply  solving  the  laminar 


boundary  layer  equations^, 71.  Results  of  these  methods  have  been  found  to  he  good  althouqh  the  prediction 
of  laminar  separation  is  not  always  as  good  as  one  would  like.  Once  the  laminar  houndary  layer  separates, 
the  calculations  become  difficult  because  the  transition  of  this  unstable  shear  layer  to  turbulent  flow 
proceeds  rapidly.  The  boundary  layer  equations  can  be  used  to  calculate  flows  with  thfn  separation  regions. 
This  procedure  is  usually  referred  to  as  the  "viscid-inviscid"  interaction  All  analytical  and 
numerical  techniques  presently  available  suffer  from  our  Incomplete  understanding  of  the  transition  process. 
Although  transition  in  an  attached  boundary  layer  can  occur  in  this  Reynolds  number  ranged.*.  at  chord 
Reynolds  number  above  about  300,000),  transition  after  separation  is  more  common. 

2.3  Separation  Bubble  and  Transition 

The  leading  edge  separation  bubble  as  sketched  in  Figure  12  is  formed  when  the  laminar  boundary  layer 
separates  from  the  surface  as  a  result  of  the  strong  adverse  pressure  gradient  downstream  of  the  point  of 
minimum  pressure,  shown  in  Figure  13.  This  separated  shear  layer  is  very  unstable  and  transition  usually 
begins  (point  T‘  on  Figure  13)  a  short  distance  downstream  of  separation,  as  a  result  of  the  amplification 
of  velocity  disturbances  present  immediately  after  separation.  After  complete  transition  from  laminar  to 
turbulent  flow  (point  T  on  Figure  13),  the  large  turbulent  shear  stresses  energize  the  shear  layer  hy 
entraining  fluid  from  the  external  stream  so  that  it  grows  rapidly  toward  the  airfoil  surface,  causing  the 
pressure  to  rise.  Reattachment  occurs  when  the  pressure  is  nearly  equal  tn  the  value  for  the  turbulent, 
boundary  layer  over  the  airfoil  with  no  separation  bubble  present,  as  shown  in  Figure  13.  The  inviscid  flow 
solution  value  of  pressure  is  also  frequently  used  to  determine  reattachment  location  since  it  approximates 
the  turbulent  boundary  layer  case.  The  region  between  separation  and  reattachment  is  referred  to  as  the 
separation  bubble.  The  fluid  in  the  laminar  portion  of  the  bubble  moves  very  slowly,  while  the  fluid  in  the 
turbulent  portion  moves  vigorously  in  a  recirculating  pattern.  Those  factors  which  affect  boundary  layer 
separation  also  affect  the  separation  bubble  and  transition  in  the  separated  shear  layer,  namely:  thickness 
of  the  boundary  layer  at  separation,  angle  of  attack,  free  stream  turhulence  level  and/or  other  free  stream 
disturbances  and  surface  roughness. 

The  transition  process  in  the  separated  shear  layer  is  the  keystone  which  determines  the  size  and  shape 
of  the  bubble,  as  well  as  how  rapidly  the  developing  turbulent  boundary  layer  grows  over  the  remaining 
portion  of  the  airfoil.  Although  the  transition  process  in  attached  shear  layers  (especially  on  flat 
plates)  has  received  a  great  deal  of  attention  (e.g.,  References  38  and  77-81)  transition  In  separated  shear 
layers  (especially  in  airfoil  type  pressure  gradients)  has  received  much  less  attention.  Almost  all  of  the 
previous  research  related  to  separated  shear  layer  transition  has  been  performed  usinq  either  a 
two-diV-MSional  backstep  geometry  82-85  or  the  exit  of  an  axisymmetric  free  jet8fi-8R.  Roth  of  these 
geometr'es  produce  a  separated  shear  layer  by  an  abrupt  disappearance  of  the  wall  over  which  the  shear  layer 
is  growing;  that  is,  separation  is  forced  hy  the  abrupt  change  in  geometry.  In  the  case  of  the  hackstep, 
reattachment  takes  place  downstream  and  a  separation  bubble  is  formed.  For  the  axisymmetric  jet 
experiments,  no  surface  or  wall  is  availahle  for  reattachment  and  no  separation  bubble  is  formed.  The 
leading  edge  separation  bubble  on  an  airfoil  is  quite  d-fferent  from  hnth  of  these  sincp  separation  is 
caused  by  an  adverse  pressure  gradient.  This  separation  bubble  and  the  separated  shear  layer  transition 
appear  to  be  more  dependent  on  the  interaction  between  the  airfoil  boundary  layer  and  the  external  flow 
field  around  the  entire  airfoil.  It  is  clear  that  while  there  must  he  some  similarities  in  the  transition 
process,  no  matter  how  it  is  caused,  there  must  also  he  significant  differences*^. 

It  is  generally  agreed  that  transition  from  laminar  to  turbulent  flow  may  he  described  as  a  series  of 
events  which  take  place  more  or  less  continuously,  depending  on  thp  flow  problem  studied.  Since  turbulence 
is  essentially  a  three-dimensional  phenomenon,  the  breakdown  of  a  two-dimensional  laminar  flow  may  hp  viewpd 
as  the  process  whereby  finite  amplitude  velocity  fluctuations,  or  traveling  wavp  disturbances,  acquirp 
significant  three-dimensional  i  ty^-vl .  The  velocity  fluctuation  or  traveling  wave  front  which  is  Initially 
straight  develops  spanwise  undulations  that  are  enhanced  by  second  order  effects,  as  depictpd  in  Flqure  14, 
Transition  has  been  very  graphically  described  as  the  process  hy  which  the  straiqht  and  parallel  vortex 
lines  of  a  two-dimensional  laminar  flow  deform  into  a  constantly  changing  and  twisting  three-dimensional 
mess  call  "turbulence"'”  .  This  Is  best  described  by  a  quote  from  Reference  78: 

"It  is  not  the  mere  presence  of  vorticity  that  characterizes  turhulence.  It  's  the  complexity 
of  the  vorticity  field.  In  a  laminar  boundary  layer,  the  vortex  lines  are  parallel  and  stacked 
near  the  wall,  like  uncooked  spaghetti.  In  the  turbulent  layer,  the  vortex  lines  are  constantly 
changing  and  twisting.  Near  the  wall,  major  entanglements  appear,  and  the  vortex  lines  may 
develop  knots  and  crossover  points.  The  spaghetti  is  cooked." 


Still  photographs  and  high  speed  movies  of  the  smoke  filaments  (i.*,*.,  streaklines)  produced  by  passing 
an  electrical  current  through  a  5  mil  wire  coated  with  oil  (Reference  R?),  clearly  delineate  the  fine 
details  of  the  separated  shear  layer  from  the  leading  edge  of  the  airfoil.  Figure  15  was  obtained  using 
this  technique  known  as  the  "smoke-wire"  method.  The  transition  pro  ess  in  the  separated  shear  layer  is 
seen  in  all  of  its  complexity.  A  preliminary  examination  of  these  smoke  photographs  substantiates  the 
notions  of  a  highly  unstable  two-dimensional  flow  which  breaks  down  in  a  very  definite  manner  to  a 
three-dimensional  turbulent  flow.  These  smoke  photographs  represent  the  most  definitive  visual  description 
of  separated  shear  layer  transition  available.  Some  structure  is  also  visiblp  in  the  devplopinq  turbulent 
flow.  Although  this  visual  technique  Is  only  suited  for  chord  Reynolds  numbers  less  than  inn,nnn,  the  basic 
transition  process  should  follow  the  same  series  of  events  at  higher  Reynolds  numbers.  For  examplp,  for  thp 
same  airfoil,  the  beginning  of  the  transition  process  moves  toward;  the  separation  location  as  the  free 
stream  velocity  is  increased.  The  length  of  the  transition  region  also  decreased  with  higher  free  stream 
velocities. 


2.4  Low  Reynolds  Number  Turbulent  Flow 


Integral  methods  o'  solving  the  turbulent  boundary  layer  equations  have  been  in  wide  use  following  the 
development  of  Prandtl's  mixing  length  hypothesis  in  192536.  Roth  the  variety  and  complexity  of  the 
integral  methods  increased  with  time  until  the  1968  Stanford  Conference  which  was  designed  to  determine  the 
accuracy  of  the  available  turbulent  boundary-layer  prediction  methods'^.  The  results  of  this  conference 
demonstrated  that  the  partial  differential  equation  methods  provided  more  accurate  predictions  of  turbulent 
boundary-layers  than  the  best  integral  techniques.  These  more  complex  methods  using  large-fast  digital 
computers  were  able  to  predict  detailed  features  of  turbulent  flows.  It  became  clear  at  this  time  that  more 
physical  experiments  were  needed  in  turbulent  flows  so  that  the  numerical  semi-empirical  turbulence  models 
could  be  improved.  These  physical  experiments  greatly  increased  our  understanding  of  f^o  basic  structure  of 

turbulence^»9!>,96  as  w>?]i  as  path  from  iamjnar  t0  turbulent  flowQ3,9R. 

The  development  of  numerical  techniques  used  to  solve  the  Navier-Stokes  equations  with  turbulence 
modelling  has  been  very  rapid  in  the  past  few  years™"!1’* .  Approximate  calculations  of  this  type,  should 
eventually  be  useful  for  engineering  purposes. 

At  low  chord  Reynolds  numbers  on  an  airfoil,  the  turbulent  boundary  layers  are  usually  in  the  staqe  of 
developing  after  attached  transition  or  redeveloping  after  reattachment  when  a  separation  huhhle  is  present. 
These  boundary  layers  are  in  general  far  from  the  fully-developed  concept.  In  these  developing  or 
redeveloping  boundary  layers  it  appears  plausible  that  their  structure  is  dependent  upon  the  mechanism  of 
transition  as  well  as  the  pressure  gradient,  etc.  Although  low  Reynolds  number  turbulent  boundary  layers  on 
airfoils  have  only  recently  begun  to  receive  a  lot  of  attention,  the  smoke  visualization  photograph  shown  in 
Figure  16  presents  a  global  view  of  this  prohlem.  Figure  16  shows  a  short  leading  edge  separation  bubble 
followed  by  a  redeveloping  turbulent  boundary  layer  on  a  NAPA  2301?  airfoil  at  a  chord  Reynolds  numhpr  of 
123,800  and  an  angle  of  attack  of  14°!^.  The  separation  bubble  (about  12*  of  the  chord  in  length)  and  the 
vortex  formation  in  the  free  shear  layer,  indicating  the  onset  of  transition,  are  clearly  visible  in  Figure 
16. 

3.  DESIGN  AND  MODELING  PROCEDURES 

It  is  not  suprising  that  the  early  airplane  designers  looked  to  flying  creatures  in  nature  for  help. 

With  birds  for  inspiration  and  the  empirical  process  of  cut  and  try  flight  and/or  wind  tunnel  experiments, 

airfoil/wing  designers  made  slow  but  steady  progress  from  the  beginning  of  this  century  until  the  early 
1920‘s.  This  procedure  is  referred  to  as  "design  by  experiment"  and  is  illustrated  in  the  left  side  of 

Figure  17.  The  additional  help  obtained  from  potential  flow  theory  in  the  1920‘s  led  to  more  systematic 

experimental  programs,  mainly  at  NACA  and  Goettingen,  which  resulted  in  catalogs  of  airfoil  sections.  These 
catalogs  of  airfoil  section  characteristics  were  very  useful  to  the  designer  183-105.  This  procedure  is 

essentially  the  same  "design  by  experiment"  one  except  trial  and  error  has  been  he’ped  a  qreaL  deal  hy  the 

use  of  potential  flow  theory. 

Although  boundary-layer  theory  was  born  at  Gottingen  in  1904,  it  wasn’t  until  the  middip  of  the  logo's 
that  it  was  developed  enough  to  he  incorporated  into  the  design  procedures  of  airfoil  sections.  Thp 
realization  of  the  importance  of  viscous  effects  resulted  in  the  well  known  NACA  6-series  "laminar  flow" 
airfoils.  The  airfoil  design  procedure,  though  guided  by  this  type  of  analysis,  was  still  mostly 

experimental  as  indicated  on  the  left  side  of  Figure  17.  The  6-series  of  airfoils  with  minor  modification 

dominated  the  aircraft  design  field  until  the  late  1950‘s  when  the  digital  computer  revolution  began  to  be  a 
factor.  At  this  point  In  history,  theoretical  considerations  could  be  used  to  do  more  than  just  guide  thp 
experimental  program. 

In  the  late  1950's,  the  power  of  large  computers  was  applied  to  designing  low  speed  airfoils  hy 
combining  potential  and  viscous  flow  theories.  In  these  so-called  "inverse  methods",  the  dpsiqner  starts 
with  the  desired  performance  characteristics,  etc.  and  ends  with  the  airfoil  shape  which  satisfies  these 
characteristics.  This  was  the  beginning  of  tailoring  airfoil  sections  for  specific  applications  and  the 
general  procedure,  referred  to  as  "design  hy  synthesis",  is  shown  on  the  right  side  of  Figure  17.  The  early 
application  of  inverse  airfoil  design  techniques  for  low-speeds  was  accomplished  primarily  by 
Wortmann38,106,107  an(j  Eppl er!08,109.  Additional  contributions  In  this  area  have  been  provided  hy  MileyllG, 
Henderson!!!*!!^,  Liebeck*13,  van  ingen114,115i  Lissamanll6,  Eppler  and  Somersll7"llq,  and  many  others. 

Although  the  particular  methods  used  to  solve  the  complex  governing  equations  for  the  inviscid  and 
viscous  flow  over  an  airfoil  in  the  inverse  mode  differ  considerably,  all  successful  inverse  design  methods 
must  contain  the  following  elements: 

Potential  Flow  Calculation  Procedure 
Laminar  Boundary  Layer  Calculation  Procedure 
Laminar  Separation  Criterion 
Laminar  Separation  Rubble  Prediction 
Transition  Criterion 

Turbulent  Roundary  Layer  Calculation  Procedure 
Turbulent  Separation  Criterion 

The  "design  by  synthesis"  approach  begins  with  the  boundary  layer  characteristics  and  their  effect  on  the 
pressure  distribution  so  that  the  aWoil  shape  which  results  meets  the  originally  desired  performance 
characteristics.  Although  several  iterations  In  the  procedure  are  usually  made  the  airfoil  shape  appears  as 
the  final  product.  This  is  exactly  the  opposite  of  the  "design  hy  experiment"  method  where  an  airfoil  shape 
is  chosen  at  the  beginning  of  the  procedure. 
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3.1  The  Eppler  Airfoil  Design  and  Analysis  Program 

The  airfoil  design  program  developed  by  Eppler  over  the  past  ?5  years  has  hepn  very  useful  in  designing 
airfoils  for  incompressible  flow  and  is  a  good  example  of  the  "design  by  synthesis"  approach.  This  program 
combines  a  conformal-mapping  method  for  the  design  of  airfoils  with  p-escrihed  velocity-distribution 
characteristics,  a  panel  method  for  the  analysis  of  the  potential  flow  about  a  qiven  airfoil,  and  a  boundary 
layer  method.  It  has  been  successfully  applied  at  chord  Reynolds  numbers  from  ?  x  10^  to  1  x  1D^. 

3.1.1  Theory 

The  airfoil  design  method  is  based  on  conformal  mapping.  This  method  differs  from  other  inverse 
methods  in  that  the  velocity  distribution  is  not  specified  for  only  one  angle  of  attack.  Instead,  angles  of 
attack  which  will  result  in  constant  velocity  over  specified  segments  of  the  airfoil  are  input.  In  other 
words,  pairs  of  parameters  are  specified:  the  segment  of  the  airfoil  and  the  angle  of  attack  relative  to 
the  zero-lift  line  which  will  result  in  constant  velocity  over  that  segment.  flf  course,  some  matching 
conditions  must  be  met  to  guarantee  a  smooth  velocity  distribution  for  all  angles  of  attack.  Toward  the 
trailing  edge,  on  both  surfaces,  a  main  pressure  recovery  can  he  specified.  Finally,  a  short  closure 
contribution  must  be  introduced  to  insure  that  the  trailing  edge  will  he  closed. 

In  reality,  the  segments  corresponding  to  the  various  input  angles  of  attack  are  not  specified  in  the 
airfoil  plane  but  rather  in  the  conformal-mapping  plane  in  which  the  airfoil  is  represented  hy  a  circle.  So 
far,  no  difficulties  have  arisen  in  correlating  the  arcs  of  the  circle  with  the  segments  of  the  airfoil. 

It  should  be  remembered  that  for  any  given  velocity  distribution  there  does  not  necessarily  exist  a 
"normal"  airfoil.  For  example,  the  closure  contributions  could  be  quite  large  which  would  result  in  a  very 
large  trai 1 ing-edge  angle.  The  closure  contributions  could  also  give  rise  to  a  region  of  negative  thickness 
near  the  trailing  edge.  Accordingly,  several  iteration  options  have  heen  included  which  allow  the 
trailing-edge  angle  to  be  specified  while  certain  input  angles  of  attack  or  the  total  amount  of  pressure 
recovery  is  iterated. 

The  potential  flow  airfoil  analysis  method  employs  panels  in  the  conformal-mappinq  plane  with 
distributed  surface  singularities.  The  geometry  of  the  panels  is  determined  by  a  spline  fit  of  the  airfoil 
coordinates,  with  the  end  points  of  the  panels  being  the  input  airfoil  coordinates  themselves.  The 
singularities  used  are  vorticities  whose  strength  are  distributed  parahol ical ly  along  each  panel.  The  flow 
condition,  which  requires  the  inner  tangential  velocity  to  he  zero,  is  satisfied  at  each  airfoil  coordinate 
(i.e.,  at  the  end  points  of  the  panels,  not  the  midpoints).  Two  angles  of  attack,  n°  and  Qfl0,  are  analyzed. 
The  flow  for  an  arbitrary  angle  of  attack  can  he  derived  from  these  two  solutions  hy  superposition.  The 
entire  procedure  does  not  require  any  restrictions  on  the  input  point  distribution,  smoothing,  nr 
rearranging  of  the  coordinates;  only  the  original  airfoil  coordinates  are  used.  An  option  is  included  by 
which  additional  points  can  he  splined  in  between  the  original  coordinates.  This  option  allows  more  precise 
results  to  be  obtained  should  a  portion  of  the  airfoil  have  a  sparse  distribution  of  points. 

A  flap  deflection  can  he  introduced  by  geometrically  rotating  part  of  the  airfoil  about  a  flap  hinqe 
point.  The  connection  between  the  forward  portion  of  the  airfoil  and  the  flap  is  defined  hy  an  arc 
consisting  of  additional  points  which  are  generated  automatically  according  to  an  input  arc  lenqth.  In 
addition,  an  option  is  included  which  allows  the  analysis  of  chord-increasing  flaps.  It  should  be  noted 
that,  while  the  airfoil  shape  which  results  form  the  exercise  of  this  option  does  have  an  increased  chord, 
it  does  not  contain  a  slot,  and  thus,  is  still  a  single-element  as  opposed  to  a  multi-element  airfoil. 

The  laminar  and  turbulent  boundary-layer  development  is  computed  using  integral  momentum  and  energy 
equations.  The  approximate  solutions  obtained  from  the  laminar  houndary-layer  method  aqree  very  well  with 
exact  solutions.  The  turbulent  boundary-layer  method  is  based  on  the  best  available  empirical 
skin-friction,  dissipation,  and  shape-factor  laws.  No  further  errors  are  introduced  by  mathematical 
simplifications  like  integrating  the  ordinary  di  "erential  equations  from  the  momentum  and  energy  laws  by 
averaging  the  right  sides  of  the  equation. 

Of  special  interest  are  the  predictions  of  separation  and  transition.  The  prediction  of  separation  Is 
determined  solely  by  the  shape  factor  based  on  energy  and  momentum  thicknesses.  (Note  that  this  shape 
factor  has  the  opposite  tendency  of  the  shape  factor  based  on  displacement  and  momentum  thicknesses.!  For  a 
laminar  boundary  layer,  there  exists  a  constant  and  reliahle  lower  limit  of  this  shape  factor,  which  equals 
1.515  and  corresponds  to  laminar  separation.  For  turbulent  boundary  layers,  no  such  unique  and  reliahle 
limit  has  heen  determined.  It  can  he  stated,  however,  that  the  turbulent  boundary  layer  will  separate  if 
the  shape  factor  goes  below  1.45  and  will  not  separate  if  the  shape  factor  remains  above  1.58.  It  has  also 
been  determined  that  thicker  boundary  layers  tend  to  separate  at  lower  shape  factors.  because  the  present 
method  yields  lower  shape  factors  for  adverse  pressure  gradients  than  other  methods,  turbulent,  separation  is 
assumed  when  the  shape  factor  equals  1.46.  The  uncertainty  is  not  as  bad  as  it  first  appears  because  the- 
shape  factor  changes  rapidly  near  separation.  Nevertheless,  results  must  he  checked  carefully  with  respect 
to  turbulent  separation. 

The  prediction  of  transition  in  the  Eppler  method  is  based  on  an  empirical  criterion  which  contains  the 
Reynolds  number,  based  on  local  conditions  and  momentum  thickness,  and  the  shape  factor.  This  criterion 
predicts  that  transition  occurs  later  if  the  shape  factor  is  higher  (i.e.,  the  pressure  gradient  is  more 
favorable).  The  criterion  also  contains  a  "roughness  factor"  which  allows  various  degrees  of  surface 
roughness  and/or  free-stream  turbulence  to  he  simulated.  The  prediction  of  transition  results  in  a  switch 
from  the  laminar  skin-friction,  dissipation,  and  shape-factor  laws  to  the  turbulent  ones,  without  chanqing 
the  shape  factor  of  the  momentum  thickness.  The  program  contains  two  options  for  fixing  transition  as  well 
as  an  option  which  allows  the  analysis  of  the  effect  of  single  roughness  elements  on  a  turbulent  houndary 
1 ayer. 


The  boundary-layer  characteristics  at  the  trailing  edge  are  used  for  the  calculation  of  the 
profile-drag  coefficient  by  a  Sqjire-Young  type  formula.  The  program  generally  predicts  slightly  higher 
drag  coefficients  than  those  measured  experimentally.  However,  the  differences  between  the  predictions  and 
experimental  measurements  depend  on  the  wind  tunnel  in  which  the  experiments  were  performed  and,  therefore, 
neither  the  skin-friction  laws  mu  any  oilier  pari  of  the  program  has  been  changed  in  orner  to  “tune"  the 
predictions  to  a  particular  wind  tunnel  or  set  of  experimental  data. 

The  lift  and  pitching-moment  coefficients  are  determined  from  the  potential  flow.  Viscous  corrections 
are  then  applied  to  these  coefficients.  The  lift-curve  slope  where  no  sepa'ation  is  present  is  reduced  to 
2 it  from  its  ti'eoretical  value.  In  other  words,  the  potential-flow  thickness  effects  are  assumed  to  he 
offset  by  the  boundary-layer  displacement  effects.  A  lift  coefficient  correction  due  to  separation  is  also 
included.  As  an  option,  the  displacement  effect  on  the  velocity  distributions  and  the  lift  and 
pitching-moment  coefficients  can  he  computed. 

3.1.2  Examples  Using  Eppler’s  Method 

The  first  example  is  a  low-Reynolds  number  airfoil  (Eppler  387)  designed  for  model  airplanes.  The 
comparison  of  the  predictions  and  the  experimental  results  of  Reference  120  is  shown  in  Figure  18.  For  a 
Reynolds  number  of  200,000  (Figure  10a),  the  agreement  between  theoretical  and  experimental  section 
characteristics  is  good.  No  pitching-moment  coefficients  were  measured.  For  a  Reynolds  number  of  10(1,000 
(Figure  18b),  the  agreement  between  theoretical  and  experimental  lift  curves  is  good.  The  predicted  maximum 
lift  coefficient  is  conservative  (low).  The  agreement  between  theoretical  and  experimental  riraq 
coefficients  is  reasonably  good  at  the  limits  of  the  low-drag  range  but  rather  poor  in  between  these  limits. 
This  poor  agreement  is  due  to  the  effect  of  laminar  separation  hubbies  or  the  measured  drag  coefficients. 
While  the  program  predicts  the  existence  of  significant  laminar  separation  bubbles  at  these  lift 
coefficients,  it  does  not  account  for  the  influence  of  the  buhhles  on  the  drag.  For  a  Reynolds  numbe1-  of 
60,000  (Figure  18c),  the  agreement  between  theoretical  and  experimental  lift  curves  is  not  as  good  as  at  the 
higher  Reynolds  numbers.  The  predicted  maximum  lift  coefficient  is  again  conservative.  The  program 
predicts  laminar  "stall"  at  moderate  lift  coefficients  (*  0.6)  and  the  turhulent  reattachment  which  results 
in  a  relatively  high  maximum  lift  coefficient.  While  the  occurrence  of  the  various  boundary-layer  phenomena 
is  predicted  well,  the  influence  on  the  drag  coefficients  is  not. 

The  second  example  is  a  natural-laminar-flow  airfoil  (NASA  NLF (1)0416)  designed  for  general  aviation 
applications.  The  comparison  of  the  predictions  and  the  experimental  results  of  Reference  121  for  a 
Reynolds  number  of  2  x  10®  is  shown  in  Figure  11.  The  agreement  between  theory  and  experiment  is  considered 
excellent.  The  magnitudes  of  both  the  angle  of  attack  for  zero  lift  coefficient  and  the  pitching-moment 
coefficients  are  overpredicted  because  the  displacement-iteration  option  was  not  exercised. 

The  third  example  is  a  flapped  natural-laminar-flow  airfoil  (NASA  NLF  ( 1 ) -021 5F )  designed  for 
high-performance  general  aviation  applications.  The  comparison  of  the  predictions  and  the  experimental 
results  of  Reference  122  for  a  Reynolds  number  of  9  x  10®  and  a  simple  flap  deflection  of  -10°  is  shown  in 
Figure  20.  The  agreement  between  theoretical  and  experimental  lift  curves  is  good.  The  pitching-moment 
«.;°ff icients  are  only  slightly  overpredicted,  even  though  the  displacement-iteration  option  was  not 
exerciseu,  •'“cause  the  displacement  effect  for  this  configuration  and  Reynolds  numher  is  small.  The 
predicted  d,  ag  tuc^icients  are  conservative  (high)  with  the  exception  of  those  at  the  lower  limit  of  the 
low-drag  range.  ’*  is  felt  that  the  disagreement  is  the  result  of  the  increased  turbulence  level  in  the 
wind  tunnel  at  1 1  ’  ligh  unit  Reynolds  number  (4.6  x  10®/ft)  which  causes  premature  transition. 

3.1.3  Remarks 

This  program  represents  a  mathematical  model  of  the  two-dimensional  viscous  flow  around  airfoils  and 
has  produced  very  good  results.  The  major  shortcoming  of  this  program  is  the  lack  of  a  separation  bubble 
model.  Since  it  appears  that  at  chord  Reynolds  numbers  below  5(10,000  a  laminar  separation  bubble  is  always 
present  and  may  vary  in  size  from  about  2%  to  20*  of  chord,  the  inclusion  of  the  bubble  in  the  analysis 
program  woula  improve  the  determination  of  lift  and  drag  for  some  cases. 

3.2  Viscous-Inviscid  Interaction  Methods 

In  the  Interest  of  predicting  laminar  separation  and  the  development  of  separation  hubbies  on  airfoils, 
another  approach  has  been  followed  by  several  researchers  72-75.  These  methods  are  based  on  the  assumption 
that  the  boundary  layer  equations  are  applicable  in  the  thin  regions  of  separated  flow  and  follow  one  of  two 
approaches. 

The  first  approach  involves  numerically  integrating  the  differential  form  of  the  boundary  layer 
equations.  In  the  attached  regions  of  the  boundary  layer,  this  can  he  accomplished  in  a  straight  forward 
manner  using  the  external  velocity  as  a  boundary  condition.  When  separation  occurs,  however,  the  equations 
become  singular.  It  then  becomes  necessary  to  use  an  Inverse  boundary  layer  solution.  This  usually  entails 
specifying  the  boundary  layer  displacement  thickness  and  solving  for  the  external  velocity.  Numerical 
instabilities  still  arise  in  the  separated  region  when  the  solution  marching  direction  is  opposite  to  the 
flow  direction.  This  difficulty  Is  usually  handled  by  making  the  so-called  FLARE  approximation,  i.e.,  by 
assuming  the  streamwise  convection  terms  are  zero.  navis  and  Carter  however,  describe  how,  in  their 
method,  the  usual  backward  difference  operator  was  switched  to  a  forward  difference  operator  In  regions  of 
reversed  flow. 

The  resulting  solution  to  the  boundary  layer  equations  must  be  verified  since  the  input  displacement 
thickness  is  only  a  guess.  Using  the  methods  of  Inviscld  flow  theory,  a  velocity  perturhatlon  due  to  the 
displacement  thickness  of  the  separated  flow  region  is  computed.  This  perturbation  is  then  added  to  the 
external  velocity  which  occurs  when  separation  Is  not  present  on  the  airfoil.  The  result  is  compared  to  the 
exter  lal  velocity  which  was  obtained  from  the  boundary  layer  solution.  If  necessary,  the  trial  displacement 
thickness  is  modified  and  the  process  is  repeated.  This  's  the  Iteration  interaction  step  for  the 


"viscous-inviscid"  procedure.  References  72  to  74  use  this  general  approach  for  obtaining  flow  solutions 
when  separation  bubbles  are  present  on  an  airfoil.  The  individual  methods  primari1y  differ  in  the  way  the 
velocity  perturbation  is  computed  and  in  the  handling  of  the  transitioning  and  turbulent  boundary  layer. 

The  method  of  Cebeci  and  C 1  ark 72  computes  the  velocity  perturbation  using  a  surface  singularity 
distribution  to  simulate  the  displacement  effects  of  the  boundary  layer.  It  apparently  can  account  for  the 
change  in  circulation  due  to  visrous  effects.  The  transition  location  is  determined  using  an  empirical 
relationship  which  is  based  on  the  flow  conditions  at  the  separation  point.  The  effects  of  freestream 
turbulence  do  not  seem  to  be  accounted  for.  The  turbulent  stress  in  the  momentum  equation  is  determined 
using  a  two  layer  turbulence  model.  In  the  transition  region,  the  stress  is  modified  by  an  intermi ttancy 
factor. 

The  technique  of  Kwon  and  Pletcher7^  uses  a  source  distribution  to  account  for  the  displacement 
thickness.  The  change  in  circulation  due  to  flow  separation  is  not  determined.  The  point  where  transition 
begins  in  the  free  shear  layer  is  assumed  to  coincide  with  the  end  of  the  separation  bubble's  constant 
pressure  region.  This  location  is  obtained  using  an  empirically  derived  relation  based  on  the  separation 
Reynolds  number.  Freestream  turbulence  effects  are  not  considered.  A  method  similar  to  that  of  Cebeci  and 
Clark7?  -js  used  to  obtain  the  turbulent  stress. 

The  computer  code  of  Davis  and  Carter7**  also  uses  a  source  distribution  to  compute  the  velocity 
distribution.  Viscous  effects  on  the  airfoil's  circulation  are  not  determined.  The  Reynolds  stress  is 
calculated  in  the  transition  and  turbulent  regions  using  a  methodology  which  accounts  for  the  magnitude  of 
the  freestream  turbulence.  This  method,  however,  was  unable  to  predict  transition  early  enough  at 
conditions  of  low  freestream  turbulence  for  separation  hubbies  to  form,  a  contradiction  to  experimental 
observations. 

The  output  of  these  calculations  is  usually  presented  in  the  form  of  wal'  shear  distributions ,  velocity 
profiles,  pressure  distributions,  and  displacement  thickness  plots.  The  results  are  encouraqinq,  however, 
their  present  application  to  the  case  of  low  Reynolds  number  flows  is  somewha  doubtful.  One  reason  for 
this  is  that  no  reliable  sensitivity  to  the  magnitude  of  freestream  turbulence  has  heen  huilt  into  the 
methods.  Furthermore,  the  programs  (except  that  of  Reference  72,  perhaps)  cannnt  handle  the  fairly  large 
separation  regions  that  occur  near  the  trailing  edge  of  airfoils  operating  at  low  Reynolds  numbers. 
Finally,  it  was  noticed  in  Reference  74  that  between  60  and  1(10  iterations  were  required  to  solve  the  case 
of  a  small  (about  3  percent  of  the  chord  in  length)  separation  hubhle.  At  low  Reynolds  numbers,  the  lenqth 
of  a  bubble  may  be  10  times  greater.  Thus,  computational  time  might  he  prohibitive  for  such  cases. 

The  second  numerical  approach,  which  is  described  hy  Gleyzes  et  al.75,  involves  solving  the  integral 
form  of  the  boundary  layer  equations  in  an  inverse  mode.  Again,  viscous-inviscid  interaction  is  required. 
In  the  laminar  and  turbulent  regions  of  the  boundary  layer,  closure  is  obtained  usinq  relationships  derived 
from  the  self-similarity  class  of  solutions.  In  the  transition  region,  the  characteristics  of  two 
fictitious  boundary  layers  one  laminar,  the  other  turhulent  are  combined  using  an  intermittancy  weighting 
function.  The  onset  of  transition  is  predicted  using  a  relation  between  an  amplification  parameter  derived 
from  the  neutral  stability  curves  of  separated  flows  and  the  local  turbulence  level.  The  predictions  of 
this  method  have  been  compared  with  experimental  data75  In  which  the  chord  Reynolds  numbers  ranged  from  ?  to 
8  million.  Turbulence  levels  were  generally  low  (less  than  (l.R  percent),  however,  a  case  in  which  the 
turbulence  level  was  2.5  percent  was  also  examined.  The  predictions  were  good  with  regard  to  the  bubble's 
effect  on  the  downstream  momentum  thickness.  The  bubble's  influence  on  the  velocity  distribution  was  not 
estimated  as  reliably.  Because  the  program  cannot  account  for  trailing  edge  separation  effects,  its  utility 
for  low  Reynolds  number  cases  is  doubtful. 

Little  work  appears  to  have  been  done  concerning  long  separation  hubbies.  This  is  probably  due  to  two 
reasons.  First,  their  negative  effect  on  performance  causes  the  airfoil  designer  to  avoid  developing 
sections  on  which  long  bubbles  can  form.  Thus,  it  is  of  primary  concern  that  one  he  able  tu  anticipate  the 
conditions  under  which  long  bubbles  form.  Knowledge  of  their  characteristics,  however,  is  not  important 
provided  that  their  formation  can  be  eliminated. 

The  second  reason  is  that  the  characteristics  of  long  bubbles  may  be  more  difficult  to  predict  than 
those  of  short  bubbles.  Their  presence  on  an  airfoil  significantly  alters  its  pressure  distribution  which 
invalidates  the  use  of  a  small  velocity  perturbation  to  account  for  the  effect  of  a  bubble.  In  addition, 
McCullough  and  Gault^®  reported  the  presence  of  a  pressure  gradient  normal  to  the  airfoil  surface  in  the 
vicinity  of  a  long  bubble.  Thus,  one  of  the  assumptions  of  the  boundary  layer  concept  is  violated.  For 
these  reasons,  the  numerical  methods  described  may  not  he  usable  for  estimating  the  properties  of  long 
bubbles. 

The  recent  work  of  Cheng7®  in  this  area  of  viscid-inviscid  interaction  methods  warrants  special 
mention.  Cheng  uses  the  triple-deck/Ki rchoff-wake  model  of  Sychev,  Messiter  and  Smith  and  its  extension  to 
represent  the  steady-state  Navier-Stokes  solutions  in  the  limit  as  RP  ■+  ».  Although  the  triple  deck  methods 
have  been  successful  in  solving  a  variety  of  laminar  flow  problems,  they  have  experienced  difficulties  in 
handling  transitional  and  turbulent  problems  because  of  our  incomplete  understanding  and  modelling  of  these 
phenomena.  This  recent  work  treats  symmetrical  and  asymmetrical  cases  with  lift  and  addresses  the  global 
problems  of  wake  closure,  including  laminar  reattachment  upstream  of  the  trailing  edge.  The  objective  of 
this  work  was  to  explore  the  manner  of  switching  from  a  grossly  separated  flow  to  a  fully  attached  flow  as 
the  profile  thickness  varied.  The  results  of  this  steady-state  approach  produced  multiple  solutions 
representing  flow  structures  with  open  and  closed  wakes.  The  investigation  of  flow  asymmetry  has  added  to 
the  multiplicity  of  the  bifurcating  solutions.  The  availability  of  two  solutions  for  an  airfoil  at 
incidence  may  explain  the  existence  of  the  upper  and  lower  boundaries  for  lift  hysteresis.  Furthermore,  the 
corresponding  result  for  a  symmetric  airfoil  at  zero  incidence  provides  an  explanation  for  asymmetric  flow 
over  geometrically  symmetrical  flow  structures  that  have  been  well  documented  in  low  Reynolds  number  airfoil 
experiments.  Although  transition  to  turbulent  flow  is  an  important  factor  in  the  wind  tunnpl  experiments, 
the  work  of  Cheng  indicates  a  completely  laminar  mechanism  for  airfoils  in  the  "suhcritical  range".  This 


theoretical  approach  (limited  though  it  may  he)  appears  to  he  the  only  me  which  predicts  the  existence  of 
hysteresis  on  an  airfoil  at  incidence  and  the  possibility  of  flow  asymmetries  for  a  symmetrical  airfnil  at 
zero  incidence.  Although  they  are  not  design  methods,  the  viscous-inviscid  interaction  methods  are  helpful 
in  analyzing  specific  airfoil  problems  and  have  added  to  our  understanding  the  role  of  separated  regions  in 
airfoil  flow  fields. 

3.3  Numerical  Modeling 

As  has  been  emphasized,  the  low  Reynolds  number  flow  field  over  an  airfoil  is  frequently  dominated  hy  a 
large  laminar  separation  bubble  and  subsequent  transition  from  laminar  to  turbulent,  or  nearly  turbulent 
conditions.  The  extent  and  nature  of  separated  flow  over  the  airfoil  is  governed  hy  the  aerodynamic 
characteristics  of  this  separation  bubble  and  transitional  flow.  Hence,  any  computational  effort  to  model 
this  flow  should  take  into  account  tne  realistic  fluid  physics,  both  in  regard  to  the  separation  and 
transitional  effects.  Boundary  layer  solutions,  no  matter  how  high  the  order,  may  overlook  some  of  the 
important  physics  of  the  separated  flow.  Therefore,  efforts  have  been  underway  to  approximately  solve  the 
complete  two-dimensional  Navier-Stokes  equations  for  low  Reynolds  number  flow  over  airfoils.  In  this 
approach,  no  modeling  is  required  for  the  separated  flow  regions.  Of  course,  with  regard  to  the 
transitional  nature  of  the  flow,  any  calculation,  Reynolds  averaged  Navier-Stokes  or  otherwise,  is  going  to 
be  dependent  on  the  particular  turbulence  model  employed.  Although  these  approximate  numerical  procedures 
are  not  suitable  at  the  present  time  for  design  purposes,  they  are  becoming  increasingly  useful  for 
analysis.  As  the  turbulence  modeling  improves  and  computational  time  decreases,  these  methods  will  be  more 
and  more  helpful  to  the  designer.  An  example  of  this  approach  is  the  work  of  Kothari,  Anderson  and 
Raghaven69.  Their  current  calculations  treat  compressible  flow  to  allow  application  from  low  to  transonic 
flight  conditions. 

3.3.1  Physical  Problem  and  Governing  Equations 

Consider  a  two-dimensional,  compressible  viscous  flow.  The  Reynolds-averaged,  time-dependent 
Navier-Stokes  equations-^  in  cartesian  coordinates  can  he  written  in  non-dimensional  formic  as 
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In  the  above,  all  quantities  are  nondimensional ized  by  reference  quantities  (subscript  r)  and  are  related  to 
the  dimensional  quantities  (primed  variables)  as 
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where 

E'  =  total  internal  energy  per  unit  mass  and 
H'  =  total  enthalpy  per  unit  mass  =  E'  +  RT'. 
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These  equations  are  transformed  Into  curvilinear  coordinates  (  5,  n)  which  are  related  to  (x,y)  through 


5  =  e(x,y)  and  n  =  n(x,y). 
The  transformed  equations  are 
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The  curvilinear  coordinates  allow  a  boundary-fitted  coordinate  system  to  he  wrapped  around  the  airfoil.  The 
details  of  the  transformation  and  the  resulting  equations  are  given  in  References  fib  and  123. 

The  molecular  viscous  coefficient  is  given  by  Sutherland's  law 


Tr  +  110 

{v  +  110  ] 


and  the  thermal  conductivity  is-  obtained  from  k'  =  u'Cp'/Pr  assuming  a  constant  Prandtl  number  of  1.0. 
These  molecular  values  are  added  to  the  turbulent  transport  properties  obtained  with  the  Raldwin-Lomax 
turbulence  model  124  which  only  includes  the  Reynolds  stress. 

The  grid  utilized  for  the  calculations  presented  here  is  generated  by  the  elliptic  grid  technique  of 
Thompson  125,  for  subsonic  applications,  this  grid  is  extended  a  large  distance  away  from  the  body,  as 
shown  in  Figure  21.  Here,  the  airfoil  is  discernible  as  a  small  spec.:  in  the  middle  of  the  grid.  A  detail 
of  the  grid  In  the  near  vicinity  of  the  airfoil  is  shown  in  Figure  22.  The  airfoil  shown  in  Figures  21  and 
22  is  the  Hiley  airfoil  investigated  in  Reference  40. 

3.3.2  Numerical  Technique 

The  discretized  Navier-Stokes  equations  are  Solved  hy  means  of  an  implicit  time-dependent 
finite-difference  technique  patterned  after  MacCormack*2".  Reference  126  treats  the  case  of  rectangular 
coordinates.  The  modifications  necessary  for  a  curvilinear  coordinate  system  are  detailed  in  Reference  6b. 
The  time-dependent  technique  starts  with  assumed  Initial  conditions  throughout  the  flowfield  (at  time  =  01, 
and  calculates  the  flow  in  steps  of  time.  At  large  times,  the  steady-state  flow  is  approached.  It  is  this 
steady  flow  that  is  of  interest;  the  time-dependent  technique  is  simply  a  means  to  that  end. 
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3.3.3  Typical  Results 

Consider  the  laminar  flow  over  a  Miley  airfoil  at  H  =  0.5  and  Rc  =  100,000  (hased  on  chord).  The 
calculdLeu  streamlines  fur  liiis  flow  die  shown  in  Figure  23  for  zero  degrees  angle-of-attack .  This 
streamline  pattern  agrees  in  general  with  smoke  flow  photographs  at  the  University  of  Notre  name  Note 
the  massive  separation  downstream  of  the  maximum  thickness  of  the  airfoil.  The  calculated  variation  of 
pressure  coefficient  over  the  airfoil  is  given  in  Figure  24,  and  is  compared  with  experimental  values  from 
Pohlen  and  Hueller^O.  Reasonable  agreement  is  obtained.  Finally,  the  calculated  variation  of  lift  and  drag 
coefficients  with  angle  of  attack  is  shown  in  Figure  25,  along  with  the  experimental  values  given  in 
Reference  40.  Again,  reasonable  agreement  is  obtained  for  a<  10°.  This  agreement  is  particularly 
important  in  light  of  the  tenuous  nature  of  low  Reynolds  number  data  for  airfoils.  Indeed,  Figure  25 
demonstrates  the  value  of  such  Navier-Stokes  calculations  for  Reynolds  numbers  low  enough  where  experimental 
data  is  questionable  or  non-existent.  The  lack  of  agreement  at  large  a  is  due  to  this  particular  set  of 
calculations  being  completely  laminar,  whereas  transition  is  occuring  in  the  actual  experiments. 

In  addition  to  the  Miley  airfoil,  another  airfoil  of  interest  viz.  a  Uortmann  FX  63-137  airfoil  was 
also  included  in  the  computational  simulation  experiment.  The  streamline  pattern  for  zero  degrees  angle  of 
attack  is  shown  in  Figure  25  for  Rc  **  100,(100.  The  flow  separates  on  both  upper  and  lower  s’des  of  thc> 
airfoil,  with  upper  surface  separation  being  more  massive.  It  should  again  be  noted  that  the  solution  here 
is  completely  laminar  and  in  the  case  of  the  actual  experiments  at  Rc  =  100,000  transition  could  occur 
causing  the  flow  on  the  lower  and  upper  surfaces  to  reattach. 

Although  approximate  numerical  solutions  of  this  type  are  not  suitable  for  design  purposes  they  make  an 
important  contribution  to  our  understanding  of  the  complex  unsteady  separated  flow  problems  encountered  with 
low  Reynolds  number  airfoils.  As  our  ability  to  model  transition  and  turbulence  improves, this  approach  will 
become  more  helpful  to  the  designer. 

4.  EXPERIMENTAL  EVALUATION  OF  DESIGN 

Although  now  airfoil  sections  may  now  be  obtained  using  the  "design  by  synthesis"  methods,  the  final 
proof  of  any  design  is  its  performance  under  actual  operating  conditions.  The  most  economical  method  of 
verifying  a  given  airfoil  design  is  to  simulate  actual  operating  conditions  (i.e.,  flight  conditions)  in  a 
wind  tunnel.  Wind  tunnel  experiments,  if  carefully  performed  and  documented,  usually  provide  a  reasonable 
evaluation  of  how  an  airfoil  section  will  perform  in  flight. 

Since  the  overall  performance  of  airfoils  at  low  Reynolds  numbers  is  very  sensitive  to  the  location  of 
transition,  care  must  be  exercised  to  understand  the  peculiar  wind  tunnel  disturbance  environment  and  its 
effect  on  transition.  Recause  each  wind  tunnel  facility  is  different  in  its  disturbance  environment,  it  may 
be  desirable  to  use  a  sailplane  or  an  appropriately  designed  RPV  as  a  test  bed  to  verify  airfoil  designs 
under  realistic  flight  conditions. 

4.1  Wind  Tunnel  Experiments 

To  evaluate  and  improve  existing  airfoil  design  procedures,  accurate  wind  tunnel  data  are  needed. 
These  data  include  lift,  drag  and  moment  measurements,  as  well  as  the  determination  of  the  location  of 
transition  and  separation  on  two-dimensional  airfoil  sections  and  finite  winqs.  The  experimental  studies 
discussed  below  will  emphasize  the  research  performed  at  the  University  of  Notre  Rarne.  The  data  presented 
are  representative  of  the  work  initiated  in  the  early  lR70’s  and  Indicative  of  the  onquing  nature  of  the 
types  of  investigations  being  carried  out  at  Notre  Dame.  As  is  typical  for  any  research  laboratory,  these 
facilities  have  grown  and  progressed  in  accordance  with  time,  resources,  and  experience. 

4.1.1  Influence  of  Experimental  Technique  and  Procedure 

The  lift  force  can  be  determined  with  acceptable  accuracy  using  a  strain  gauge  balance.  If  small 
pressure  differences  can  be  measured  accurately,  then  a  reasonably  accurate  lift  force  can  also  he  nhtained 
by  integrating  the  static  pressure  distribution  around  the  airfoil.  In  an  atmospheric  wind  tunnel  at  low 
Reynolds  numbers  this  requires  accurate  pressure  difference  measurements  below  1  mm  of  water. 

Since  airfoil  draq  forces  are  at  least  an  order  of  magnitude  smaller  than  the  lift  forces,  they  are 
much  more  difficult  to  measure  accurately  when  tests  are  conducted  at  low  Rc*  "‘ensure  these  very  small 
drag  forces  {e.g.,  0.015  N  for  a  250  nw  chord,  403  ms  span  NACA  663-018  airfoil  section  at  Rc  =  40,000  or 
0.01N  for  a  152  mm  chord,  412  mn  span  FX  63-137  airfoil  section  at  Rc  s  00,000),  a  very  sensitive,  specially 
designed  strain  gauge  balance  arrangement  with  high  signal  to  noise  ratio  and  low  drift  electronics  is 
needed.  A  balance  system  of  this  type  is  within  the  state-of-the-art  and  has  bpen  under  development  at  the 
University  «. f  Notre  Dame.  The  major  difficulty  with  this  technique  is  the  determination  of  Interference 
effects  between  airfoil  and  side  plates1?7.  A  sketch  of  some  of  the  possible  three-dimensional  flow 
phenomena  which  may  occur  In  the  vicinity  of  the  alrfoil/side  plate  region  is  shown  in  Figure  27  .  These 
phenomena  are  dominated  by  flow  separation  and  the  formation  of  a  corner  vortex.  It  is  cl  par  that,  thp  force 
contribution  of  the  end  of  a  finite  span  airfoil  is  different  from  the  infinite  span  airfoil  it  attempts  to 
represent.  At  low  chord  Reynolds  numbers  (i.e.,  Rj  <  1(10,000),  the  reqinn  affected  hy  side  plate 
interference  may  be  as  large  as  5t  of  the  span  ( 1  OX  of  the  span  If  sidp  plates  are  at  both  ends  of  the 
airfoil).  As  Reynolds  number  ( persists ,  this  region  nf  ihree-rl  i  ••eps  I  gnal  flow  decreases  in  size  and  effect. 
As  a  result  of  this  three-dimensional  corner  flow  region  near  -:he  ends  of  the  airfoil,  the  measured  drag 
forces  are  expected  to  be  higher  than  for  the  infinite  span  airfoil. 

Another,  often  used,  method  of  determining  drag  is  the  wake  traverse  method.  This  method  has  been  used 
successfully  at  high  Reynolds  numbers  for  several  decades.  At  first  glance  there  appears  to  be  no  reason 
why  the  wake  traverse  should  not  produce  good  results  at  low  Reynolds  numbers  if  the  static  and  total 
pressures  can  be  measured  accurately.  However,  for  Rc  <  10(1,00(1,  most  airfoil  wakes  are  composed  of 
large-scale  vortices  which  produce  an  unsteady  or  oscillating  wake  somewhat  similar  to  the  one  behind  a 
circular  cylinder  (see  Figure  28). 
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A  fixed  static  and  total  pressure  rake,  as  indicated  in  Figure  28,  is  subject  to  errors  related  to  the 
changing  flow  direction.  The  velocity  distribution  obtained  from  the  wake  traverse  is  used  in  the  momentum 
equation  written  in  the  direction  parallel  to  the  test  section  centerline.  The  airfoil  drag  is  then  assumed 
to  be  equal  to  the  decrease  in  momentum  in  this  direction.  Since  these  low  Reynolds  number  airfoil  flows 
are  dominated  uy  large-scale  vortices,  accurate  measurement  of  the  velocity  component  parallel  to  the  test 
section  centerline  is  very  difficult,  if  possible  at  all,  with  a  rake  wake.  The  following  quote  from  Pope 
and  Harperl28  supports  this  argument:  "The  wake  survey  cannot  be  used  to  measure  the  drag  of  stalled 
airfoils  or  of  airfoils  with  flaps  down."  The  drag  determined  hy  this  method  appears  to  he  lower  than  the 
actual  drag.  Furthermore,  at  these  low  Reynolds  numbers,  laminar  separation  and  transition,  usually  on  the 
airfoil  upper  surface,  very  often  produce  a  large-scale  spanwise  flow  structure^7  (see  Fiqure  PR).  The 
measurement  of  static  and  total  pressures  with  a  rake  in  such  an  oscillating  wake  can  be  subject  to 
considerable  error,  even  if  the  small  pressure  differences  can  be  accurately  resolved.  This  problem  is 
compounded  when  a  significant  spanwise  flow  structure  is  present.  To  adequately  account  for  the  spanwise 
variations,  the  rake  must  be  traversed  in  the  spanwise  direction. 

A  recent  study  of  the  variation  of  the  measured  drag  in  the  spanwise  direction  using  a  wake  rake  was 
performed  by  Althaus^9.  Using  an  integrating  rake,  an  analog-digital  converter  and  a  digital  computer 
system,  drag  data  for  1,000,000  <  Rc  <  3,000,000  was  ohtained.  The  airfoil  which  spanned  the  tunnel  had  a 
span  of  73  cm  and  a  chord  of  50  cm.  The  results  of  this  experiment  are  shown  In  Figure  30.  Drag  data  were 
obtained  over  the  center  30  cm  of  the  airfoil  span.  Although  the  lowest  chord  Reynolds  number  studied  was 
1,000,000,  the  spanwise  drag  variation  using  this  technique  is  significant.  The  spanwise  flow  structure 
increases  in  scale  as  the  Reynolds  number  decreases  and,  thus,  the  maximum  deviation  from  the  mean  in  the 
measured  drag  increases. 

In  view  of  the  large-scale  vortices  produced  in  the  wake  and  the  spanwise  flow  structure,  fhe  strain 
gauge  balance  method  appears  to  he  attractive  if  the  side  plate  interference  effect  could  he  determined  or 
el iminated. 

4. 1.1.1  Wind  Tunnel  and  Force  Ralance 

A  large  number  of  experiments  have  been  performed  in  the  non-return,  low  speed  wind  tunnels  in  the 
Aerospace  Laboratory  of  the  University  of  Notre  Flame.  These  tunnels  are  capable  of  producing  low  turbulence 
intensities  over  the  normal  range  of  tunnel  velocities,  9  to  28  m/sec  (see  Figure  31).  The  turbulence 
intensity  which  varies  only  slightly  over  this  range  has  been  found  to  he  equal  to  or  less  than  ahout  O.lt*. 
Twelve  anti-turbulence  screens  precede  the  contraction  cone  which  has  a  square  cross-section  and  a 
contraction  ratio  of  24:1.  The  test  section,  which  is  610  mm  square  (24  in  x  24  in)  and  1828  mm  (72  in) 
long,  and  diffuser  are  separated  hy  a  lfll  mm  wide  section  of  foam  rubber,  used  to  minimize  vihrations  from 
the  fan  motor.  The  eight-bladed  fan  is  driven  by  a  15  hp  A0  motor  with  a  variahle  speed  drive.  The  motor 
and  fan  assembly  is  mounted  outside  the  laboratory  in  a  protected  structure  isolated  from  the  diffuser.  The 
lower  limit  of  the  tunnel's  velocity,  normally  9  m/sec,  can  he  extended  to  as  low  as  2.1  m/sec  hy  the 
addition  of  one  or  two  flow  restrictors  between  the  vibration  insulation  and  the  test  section.  The  flow 
restrictors  are  made  from  ordinary  plastic  drinking  straws  5  mm  1. 8.  by  200  mm  long,  and  also  serve  to  damp 
out  slight  surges  in  tunnel  velocity  caused  hy  outside  wind  gusts. 

The  test  section  used  to  obtain  lift  and  drag  force  data  had  an  externally  mounted,  two-component 
strain  gauge  balance  with  a  dual  flexure  system  for  lift  and  drag.  For  low  loads  as  small  as  0.01  N  (0.04 
oz),  a  very  sensitive  flexure  was  engaged;  a  stiffer  flexure  was  engaged  at  around  12.74  N  (46  oz.).  The 
model  was  mounted  with  the  span  vertical  between  the  two  square,  smooth  plates,  9.5  x  610  x  610  mm(3/8  x  24 
x  24  in.),  bolted  to  the  inside  of  the  test  section  so  that  the  model  and  the  sting  "floated"  between  the 
plates  as  shown  in  Figure  32.  The  gaps  between  the  model  and  the  endplates  were  hpld  as  close  to  0.51  mm 
(0.020  in.)  as  possible  to  minimize  leakage  of  flow  through  them  without  risking  contact  hetwepn  model  and 
plate.  The  configuration  thus  approximated  on  airfoil  aspect  ratio  of  Infinity.  The  modpl  stinq  was 
shielded  from  the  air  flow  hy  a  streamlined  covering.  Thp  anqle  of  attack  of  the  airfoil  was  changed  with 
the  tunnel  running  hy  means  of  an  electric  motor/gear  arrangement.  Changes  in  the  angle  of  attack  could  he 
made  to  within  *  0.05  degrees. 

To  investigate  the  intei  action  hetween  the  boundary  layer  on  the  plastic  plates  and  the  ends  of  the 
airfoil  model,  the  arrangement  shown  in  Fiqure  33  was  used.  For  this  configuration,  the  airfoil  model  (s 
made  in  three  pieces  with  the  top  and  bottom  sections  attached  to  the  plates  and  thp  center  section  free  to 
float.  The  center  section  is  attached  to  the  external  balance  with  a  sting  which  passes  through  a  hole  in 
the  top  piece  of  the  model.  The  gaps  between  the  center  section  of  the  model  and  thp  end  sections  wptp  n.61 
mm  (0.020  in.).  To  change  angle  of  attack,  the  center  section  was  rotated  with  the  motor  drive  to  thp 
desired  angle,  after  which  the  top  and  bottom  sections  were  rotated  manually  and  secured  to  the  plates. 

The  signals  from  the  strain  gauge  balance  were  sent  directly  into  the  strain  gauge  amplifier.  The  lift 
and  drag  amplifier  gains  were  designed  to  be  different  due  to  the  larger  forces  resulting  from  lift;  thus, 
the  drag  gain  was  about  five  times  greater  than  the  lift  gain.  The  signals  were  subsequently  recorded  on  a 
strip  chart  recorder. 

The  NACA  663-OI8  and  5.64%  thick  EPPLER  61  airfoil  geometries  were  used  in  these  experiments,  as  shown 
in  Figure  34.  The  airfoil  models  were  made  in  the  Aerospace  Laboratory  at  the  University  of  Notre  name  by 
machining  the  models  out  of  aluminum  and  then  making  epoxy  molds.  Any  number  of  identical  models  with  a 
smooth,  plastic-like  finish  can  be  cast  using  each  mold.  Three  NACA  663-OIR  models  were  used,  each  having  a 
chord  of  250  mm  and  a  span  of  403  mm.  One  model  was  used  as  a  full  span  force  model  (Figure  32),  one  was 
cut  into  three  pieces  with  only  the  center  254  mm  span  attached  to  the  force  balance  to  check  the  side 
plate/airfoil  interaction  (Figure  33)  and  the  third  was  instrumented  with  98  static  pressure  taps. 


A  discussion  of  the  free  stream  disturbance  level  in  the  wind  tunnel  will  he  given  in  section  4.1.2, 


Two  smooth  EPPLER  61  airfoil  models  were  also  used,  each  having  a  chord  of  IPS  mm  and  a  span  of  inn  mm. 
One  of  these  models  was  milled  from  the  aluminum  model  and  was  cut  into  three  pieces  to  check  the  side 
plate/airfoil  interaction.  Only  the  center  274  mm  section  of  the  airfoil  was  attached  to  the  force  balance. 
Data  from  an  earlier  study^O  using  a  reinforced  balsa  wood  EPPLER  61  model  with  a  PSfl  mm  chord  and  a  424  mm 
span  was  used  for  comparison  purposes. 

4. 1.1. 2  Pressure  Measurement 

An  automated  pressure  measuring  system  was  developed  to  improve  the  speed  and  accuracy  of  measuring 
pressures  around  an  airfoil  section.  The  system  components  included  an  Apple  11  plus  mini-computer  with 
four  floppy  disk  drives,  a  printer,  and  a  video  monitor.  Installed  in  the  Apple  II  Plus  were  an  Interactive 
Structures  AI 13-12-bit  analog  to  digital  converter  data  acquisition  system  and  a  Mountain  Computer  A/n  +  n/A 
8-bit  input  plus  output  system.  Each  of  these  A/D  devices  could  accommodate  16  channels  of  input  and  the 
Mountain  D/A  could  accommodate  16  channels  of  output.  Two  Setra  Systems  electronic  manometers  with  analog 
output  were  used  (Models  3398,  calibrated  from  0  to  14  mm  water  TO  tj  0.55  in.  waterl,  and  339H,  calibrated 
from  0  to  140  inn  water  [0  to  5.5  in.  waterl)  to  measure  tunnel  dynamic  pressure  from  a  pitot  static  tube, 
and  to  measure  the  pressure  differences  at  each  airfoil  pressure  port.  Roth  electronic  manometers  used 
differential  pressure  transducers  and  had  digital  read-out  to  four  significant  figures.  These  devices  had  a 
listed  accura cy  of  0.14%  of  full  scale. 

The  pressure  ports  in  the-  model  surface  were  connected  with  tubing  to  junction  switch  wafers  in  three 
groups  of  33  ports.  The  output  of  the  wafers  was  connected  to  a  Scanivalve  body.  A  dummy  transducer  was 
installed  in  the  body  and  the  output  connected  to  one  of  the  electronic  manometers.  The  standard  solenoid 
controller  was  used  for  the  Scanivalve,  but  the  Mountain  D/A  output  was  used  in  conjunction  with  a 
transistor  trigger  circuit  to  switch  and  home  the  Scanivalve.  The  analoq  outputs  of  the  electronic 
manometers  were  connected  to  two  channels  of  the  Interactive  Structures  AI13  unit  to  measure  the  (Pi  -P„)  at 
each  port  on  the  model  minus  the  free  stream  static  pressure  and,  simultaneously,  the  free-stream  dynamic 
pressure  (Q„). 

The  pressure  tap  model  had  ports  situated  on  the  leading  edge  and  near  the  trailinq  edqe,  with  46  ports 
on  both  surfaces  along  the  chord.  Four  extra  ports  were  placed  spanwise  on  the  model  at  one  chordwise 
location,  giving  a  total  of  five  ports  spaced  across  the  span  for  two-dimensional  flow  checks.  The  port 
tubing  was  first  connected  to  the  wafers  which  reduced  the  number  of  tubes  so  that  the  5canivalve  could  he 
used.  The  Scanivalve  was  then  connected  to  the  electronic  manometer  so  that  all  the  airfoil  ports  could  he 
easily  accessible.  The  pitot  static  tube,  mounted  ahead  of  the  upper  side  plate,  was  connected  across  the 
second  manometer;  the  static  pressure  line  was  also  connected  to  the  manometer  for  port  pressure  to  provide 
A  P. 


A  computer  program  was  written  in  the  Apple  Raslc  language.  This  program  homes  the  Scanivalve  to  its 
first  port  and  then  waits  until  the  operator  sets  the  angle  of  attack,  tunnel  dynamic  pressure  Q,  and  thp 
wafers  to  the  correct  position.  When  ready,  the  computer  takes  200  samples  of  each  manometer  voltage  over  a 
period  of  about  10  seconds,  averages  them  to  get  a  mean,  converts  the  voltages  to  inches  of  water  using  the 
manometer  calibration  constant,  and  computes  the  pressure  coefficient  Cp  *  -P^/Q  at  each  port.  After 
sampling  through  the  first  33  ports,  the  program  enters  a  wait  until  the  operator  switches  the  wafers  to  the 
next  set  of  ports.  This  procedure  continues  similarly  for  the  last  32  ports  which  completes  the  run.  The 
computed  Cp  values  and  other  pertinent  test  information  such  as  the  date  of  test,  average  Q,  angle  of 
attack,  atmospheric  pressure,  and  temperature  are  then  sent  to  the  printer  for  a  hard  copy. 

4. 1.1. 3  Discussion  of  Results 

Force  data  obtained  from  the  strain  gauge  balance  are  presented  for  the  NACA  663-OIR  and  the  FPPLFR  61 
airfoils.  Drag  force  data  obtained  by  integrating  the  static  pressure  distributions  for  the  NACA  663-OI8 
airfoil  are  also  presented.  All  force  data  are  corrected  for  solid  body  blocking  wake  blocking,  streamline 
curvature,  and  longitudinal  buoyancy  according  to  the  methods  described  by  Rogers!-^.  The  applicability  of 
these  standard  corrections  was  demonstrated  in  Reference  60  where  two  different  sized  airfoils  produced 
almost  identical  corrected  values  of  the  c£  and  Crf  at  Rc  =  87,000. 

Lift  and  drag  data  for  the  NACA  663-OI8  airfoil  were  obtained  from  the  strain  qauqe  balance  for  61,000 
<  Rc  <  300,000.  The  results  for  Rc  =  61,000  and  300,000  are  shown  ir.  Figures  35  and  36  respectively.  Data 
are  shown  on  each  of  these  fiqures  for  the  one  piece  airfoil  model  and  for  the  three  piece  airfoil  model. 
The  difference  in  the  drag  coefficient  obtained  from  the  two  models  is  attributed  to  thp  elimination  of  most 
of  the  airfoil/end  plate  interaction  for  the  three  piece  model.  A  smaller  and  less  consistent  difference  is 
noticeable  for  the  lift  coefficient  over  this  Reynolds  number  rarge.  A  comparison  of  the  minimum  draq 
coefficient  for  this  airfoil  with  the  one  piece  and  three  piece  models  is  shown  in  Fiqurp  37.  Thp 
uncertainty  in  these  data  as  obtained  by  the  method  of  Kline  and  McCl  Intock*  is  also  shown.  The 
determination  of  the  uncertainty  of  single-sample  measurements  is  particularly  important  in  this  problem 
because  of  the  difficulties  in  the  measurement  and  because  the  truedraq  force  is  not  known.  Thp  differpncp 
between  the  Cdm)n  measured  for  the  two  models  below  about  Rc  =  80,000  is  overshadowed  hy  the  uncertainty  in 
the  measurements.  As  Rc  increases  to  about  80,000  the  trend  becomes  clear.  The  difference  in  Cdmjn  at  Rc  = 
60,000  is  about  15%  and  decreases  to  about  10%  at  Rc  =  300,000.  Studies  using  a  circular  cylinder  and  this 
same  balance”^  produced  drag  coefficients  which  were  10%  higher  than  other  published  data.  Figure  37  also 
contains  four  data  points  which  were  obtained  by  integrating  static  pressure  distributions.  A  typical 
static  pressure  distribution  is  shown  in  Figure  38.  At  this  Reynolds  number  (90,000)  thp  minimum  drag 
occurred  at  an  angle  of  attack  of  7°.  Since  the  pressure  distributions  do  not  contain  end  plate  effects  and 
since  pressure  drag  is  the  predominate  source  of  drag  it  is  not  surprising  that  these  points  fall  close  to 
and  slightly  lower  than  the  three  piece  model  data. 


Uncertainty  in  force  balance  data  has  subsequently  been  reduced  hy  a  factor  of  two  by  usinq  the  Apple  II 
data  acquisition  system. 


The  EPPLER  61  airfoil  has  been  studied  extensively  at  low  Reynolds  numbers.  Lift  and  draq  data  from 
the  strain  gauge  balance  and  smoke  visualization  data  have  been  presented  in  Reference  60.  Tho  use  of  the 
standard  tunnel  corrections  was  verified  using  25  cm  and  a  12.5  cm  chord  models.  The  reproducihi 1 i ty  of 
this  type  of  data  was  verified  for  the  one  piece  model^O.  Experiments  were  performed  to  evaluate  the 
ai  rfoi 1 /endol ate  interaction  using  one  oiece  and  three  oiece  airfoil  models  in  the  chord  Revnolds  nnmhor 
range  from  36,000  to  60,000.  The  results  of  these  experiments,  including  the  uncertainties,  were  similar  to 
the  results  for  the  NACA  663-OIB  airfoils.  The  minimum  drag  coefficient  was  found  to  he  about  161  lower  for 
the  three  piece  model  as  compared  to  the  one  piece  model. 

The  EPPLER  61  data  obtained  at  the  Oelft^?  and  Stuttgart^3  Universities  are  compared  with  data 
obtained  at  Notre  Dame  in  igure  39  for  Rc  ^  30,000.  The  nelft  and  Stuttgart  data  show  slightly  higher 
Cpnax  values  and  what  appear'  to  be  an  earlier  stall  as  angle  of  attack  increases.  The  greatest  difference 
in  these  data  sets  is  in  the  (ra'  coefficient  where  the  Delft  and  Stuttgart  results  are  as  much  as  501  lower 
than  those  obtained  at  Notre  uame.  The  one  piece  airfoil  model  used  in  the  Notre  Dame  force  balance  was 
affected  by  the  endplate  boundary  layer  interaction.  Eliminating  this  ai  rfoi  1 /erdplate  effect,  as  shown 
above,  only  reduces  the  minimum  drag  coefficient  about  161.  Why  are  the  Delft  and  Stuttgart  drag,  data 
whictt  were  obtained  with  a  wake  traverse,  still  about  411  lower?  The  difficulties  in  making  wake  pressure 
measurements  in  an  unsteady  flow  with  constantly  changing  flow  direction  probably  accounts  for  some  of  this 
difference  since  the  experiments  of  Delft  and  Stuttgart  use  the  wake  rake  to  determine  draq.  In  addition, 
the  spanwise  position  of  the  wake  traverse  may  account  for  another  part  of  this  difference.  Since  the  true 
value  of  the  drag  is  not  known,  further  studies  will  be  necessary  to  completely  resolve  this  prohlem.  It  is 
clear  that  all  measurements  at  low  Reynolds  numbers  should  he  approached  with  these  prohlems  in  mind. 

4.1.2  Influence  of  Free  Stream  Disturbances 

The  disturbance  environment  present  in  the  test  section  of  a  low  speed  wind  tunnel  is  usually 
determined  by  free  stream  turbulence  (velocity  fluctuations),  acoustic  phenomena  (pressure  fluctuations)  and 
mechanical  vibrations.  The  free  stream  turbulence  level  depends  on  the  history  of  the  flow  in  the  settlinq 
chamber,  flow  straightners  or  screens  and  the  inlet  leading  to  the  test  section.  Acoustic  phenomena  are 
related  to  the  noise  emitted  from  turbulent  houndary  layers  on  the  side  walls,  unsteady  separated  flow 
regions,  and  the  fan  and  its  associated  drive  system.  Mechanical  vibrations  may  be  caused  by  riqid  couplinq 
of  the  fan  and  drive  system  to  the  wind  tunnel  as  well  as  by  the  unsteady  wakes  of  probe  and  model  supports. 
Although  these  factors  which  determine  the  disturbance  environment  may  he  reduced  and  controlled,  they 
cannot  be  completely  eliminated.  It  is  apparent  that,  in  general,  each  wind  tunnel  has  a  different 
disturbance  environment  which  is  a  function  of  its  design  and  method  of  fabrication.  Recause  airfoil 
boundary  layers  are  sensitive  to  small  disturbances,  accurate  wind  tunnel  models  are  very  important  in  the 
evaluation  of  a  given  design.  Furthermore,  because  the  forces,  pressure  differences  and  velocities  are 
small,  a  great  deal  of  care  must  be  exercised  to  obtain  accurate  and  meaningful  d-^a.  It  is  not  surprising, 
therefore,  that  similar  experiments  on  the  same  geometry  model  at  low  Reynolds  numbers  often  produce  results 
which  differ  from  one  wind  tunnel  to  the  next. 

4. 1.2.1  Experimental  Apparatus  and  Procedure 

This  research  was  conducted  in  the  Aerospace  Laboratory  of  thp  University  of  Notre  Dame.  The  equipment 
used  consisted  of  a  wind  tunnel  and  a  strain  gauge  balance  with  its  associated  electronics  for  lift  and  drag 
measurements,  electronic  manometers,  and  a  micro-computer  system  as  described  earlier.  The  uncertainty  of 
the  earlier  lift/drag  measurements  was  reduced  by  half  by  using  the  Apple  II  data  acquisition  system  instead 
of  the  strip  chart  recorder.  A  hot-wire  anemometer  and  sound  pressure  level  analyzer  were  used  to  analyze 
the  test  section  environment.  The  flow  visualization  studies  utilized  the  smoke  tube  method,  strobe  liqhts, 
and  photographic  equipment  described  by  Mueller^.  The  airfoil  sections  used  for  thesp  studips  were  thp 
Lissaman  7769  profile  described  by  Lissaman^5  and  Rurke^^  and  the  Miley  MD6-13-12R  profilers.  Roth  of 
these  airfoil  sections  were  designed  for  a  chord  Reynolds  number  of  approximately  Rnn,DDD,  the  Lissaman  for 
the  Gossamer,  Condor  and  Albatross  man-powered  aircraft  and  the  Miley  for  the  inhoard  portion  of  a 
helicopter  rotor.  These  two  airfoils  performed  according  to  their  designs  near  Rc  =  600,000.  The 
geometries  of  these  airfoil  sections  are  shown  in  Figure  40. 

The  Lissaman  7769  airfoil  model  used  in  the  experiments  was  constructed  of  wood  using  two  steel  end 
plates  machined  to  the  profile  of  the  airfoil.  The  wood  was  coated  with  an  epoxy  and  finished  to  qivp  a 
smooth  surface.  The  airfoil  model  had  a  437  mm  span  and  a  249  mm  chord.  Fxperiments  were  also  performed 
using  two  Miley  airfoil  models  cast  from  the  same  mold.  These  smooth  epoxy  models  each  had  a  chord  of  25D 
mm  and  a  span  of  421  mm.  Dne  model  was  used  as  a  force  model  as  well  as  for  flow  visualization  and  hot-wire 
experiments.  The  other  model,  cast  with  90  static  pressure  taps,  was  used  for  pressure  measurements. 

A  hot-wire  anemometer  system  was  used  to  determine  mean  tunnel  velocities  and  turbulence  intensities  in 
the  test  section  under  varying  test  conditions.  All  hot  wire  measurements  were  made  with  a  nisA  Type  55P11 
hot-wire  probe,  having  a  sensing  element  5  microns  in  diameter  and  1.5  mm  long.  The  hot-wire  anemometer 
system  consisted  of  the  DISA  55M10  constant-temperature  anemometer  which  could  be  tuned  to  respond  to 
frequencies  up  to  50,000  Hz  with  a  hot-wire  temperature  between  250°C  and  300°C.  A  DISA  55D1D  linearizer 
was  used  and  adjusted  so  that  the  voltage  output  of  the  anemometer  corresponded  directly  to  the  tunnel 
speed.  Previous  studies  conducted  at  the  Notre  Dame  Aerospace  Laboratory  by  Kegelman^?  showed  that  the 
linearizer  produced  a  small  amount  of  electronic  noise  which  became  a  significant  part  of  the  total 
turbulent  signal  at  turbulence  intensities  below  0.11.  The  output  of  the  linearizer  was  monitored  on  a  Data 
Precision  Corporation  V-45  Digital  voltmeter  to  obtain  tunnel  velocities.  Simultaneously,  the  output  signal 
was  filtered  using  a  1  Hz  high  pass  filter  and  the  remaining  AC  component  was  amplified  to  produce  a  useful 
voltage  on  a  T.S.I.  Model  1076  true  RMS  meter,  which  (knowing  the  amplification)  yielded  thp 
root-mean-square  turbulence  Intensity,  urms  x  100/11.  When  taking  data  for  frequency  analysis,  the  niSA 
linearizer  was  bypassed  in  order  to  eliminate  any  60  cycle  noise  from  the  frequency  spectrum.  All  analog 
voltage  outputs  were  sampled  using  the  A/D  capabilities  of  thp  Apple  II  micro-computer  and  associated 
software  developed  specifically  for  these  experiments. 


Sound  pressure  level  measurements  were  accomplished  through  the  use  of  a  Rruel  and  Kjaer  Frequency 
Analyzer,  type  2107.  The  type  2107  is  an  AC  operated  audio  frequency  analyzer  of  the  constant  percentage 
bandwidth  type.  Although  designed  as  a  narrow  band  sound  and  vibration  analyzer,  it  may  he  used  for  any 
kind  of  frequency  analysis  within  the  ranqe  of  20-20.000  Hz.  The  freauencv  analyzer  was  combined  with  a 
Bruel  and  Kjaer  Level  Recorder  Type  2305,  allowing  frequency  amplitude  diagrams  to  he  recorded  automatically 
on  pre-printed  frequency-calibrated  paper.  The  Selective  Amplifier  Section  of  the  sound  pressure  level 
equipment  was  used  as  a  narrow-band  analyzer  which  could  he  continuously  varied  from  20  to  20,000  Hz.  The 
octave  selectivity  was  set  at  45  decibels  to  produce  the  narrowest  possihle  bandwidth  of  approximately  1/3 
octave.  Two  iifferent  condenser  microphones  were  used  in  the  analysis.  All  of  the  testing  was  accomplished 
with  the  use  of  a  Bruel  A  Kjaer  Nose  Cone  1IA0386  designed  to  reduce  the  aerodynamical ly  induced  noise 
present  when  the  microphones  are  exposed  to  high  wind  speeds.  The  nose  cones  were  designed  to  replace  the 
normal  protection  grid  of  the  microphone  cartridge  and  were  of  a  highly  streamlined  shape  with  a  hiqhly 
polished  surface.  A  fine  wire  mesh  around  the  circumferance  of  the  cone  allowed  sound  waves  to  penetrate  to 
the  microphone  diaphragm.  As  most  of  the  tunnel  noise  was  propagated  upstream  from  the  fan  blade  the 
directional  characteristics  of  the  microphone  were  an  important  consideration.  Rruel  and  Kjaer^fi  descrihe 
the  omnidirectional  characteristics  of  the  microphone  when  the  nose  cone  is  used.  Any  body  regardless  of 
how  streamlined  produces  some  aerodynamical ly  induced  noise.  For  the  purposes  of  this  study  it  was  assumed 
that  the  aerodynamic  noise  produced  by  the  microphone  with  its  nose  cone  would  not  exceed  the  aerodynamic 
noise  produced  by  any  similarly  streamlined  body  and  the  results  should  thus  give  a  reasonable  indication  of 
the  acoustic  environment. 

All  measurements  were  made  in  the  South  Tunnel  at  the  University  of  Notre  name  Aerospace  Laboratory. 
The  force  balance  was  used  to  collect  lift  and  drag  data  on  the  Lissaman  airfoil  over  an  angle  of  attack 
range  from  -2f  °  to  25°.  The  tunnel  was  always  started  with  the  airfoil  at  0°  angle  of  attack.  All  data 
were  correcte.:  utilizing  the  AGARD  methods  described  earlier^R.  Once  the  tunnel  was  on,  the  anqle  of 
attack  was  changed  to  -20°  and  the  experiment  was  conducted  increasing  the  angle  of  attack  to  +25°  then 
decreasing  it  back  to  0°.  During  the  entire  test  the  tunnel  was  left  running.  At  each  angle  of  attack  the 
microcomputer  sampled  the  lift,  drag,  and  the  dynamic  pressure  simultaneously  over  a  5  second  period  takinq 
100  samples  to  be  averaged.  Zero  lift  and  drag  voltages  were  measured  just  prior  and  immediately  after  the 
experiment  to  account  for  any  amplifier  drift.  The  data  was  corrected  for  this  drift  based  on  an  assumed 
linear  drift  during  the  period  of  the  experiment.  Calibration  of  the  force  balance  was  checked  every  fourth 
test.  The  force/voltage  calibration  remained  constant  over  long  periods  of  time. 

Both  the  hot-wire  and  sound  level  measurements  (fiR  cm  downstream  from  the  entrance  to  the  test  section) 
were  taken  at  a  location  8  cm  ahead  of  the  leading  edge  of  the  airfoil.  For  convenience  the  model  was 
removed  when  free-stream  measurements  were  made.  Figure  41  shows  the  experimental  set-up,  includinq  the 
location  of  the  flow-restrictor  and  the  turbulence  screen. 

Previous  studies  by  Jansen^?  and  Kegelman^7  have  shown  that  free-stream  turbulence  and  sound  levels  at 
a  given  streamwise  location  are  invarient  across  the  test  section  except  at  locations  near  the  wall.  In 
this  study  the  probes  were  centered  in  the  test  section  between  the  sldp-plates  at  a  heiqht  comparable  with 
the  location  of  the  leading  edge  of  the  airfoil  model. 

Measurements  were  made  over  a  velocity  range  of  2.5  to  31  m/s  with  different  combinations  of  thp 
turbulence  screen  and  flow  restrictors.  The  hot-wire  anemometer  was  used  to  monitor  frpp-strpam  velocity, 
measure  free-stream  turbulence  Intensities,  and  determine  frequency  content  of  thp  turbulent  velocity 
fluctuations.  Frequency  spectra  were  obtained  by  sampling  the  hot-wire  anemometer  signal,  acquiring  1024 
samples  at  a  specified  sampling  rate.  The  sampling  rate  could  he  selected  to  Increase  thp  frequency 
resolution  at  low  frequencies  or  to  increase  the  overall  frequency  range  capability.  Two  frequency  spectra 
were  taken  at  each  testing  point,  one  over  a  range  of  0  to  500  Hz  with  a  resolution  of  0.9R  Hz  and  the  other 
over  a  range  of  0  to  5000  Hz  with  a  resolution  of  9.8  Hz. 

The  Bruel  and  Kjaer  frequency  analyzer  was  used  to  determine  total  sound  pressure  levels  in  the  test 
section.  Calibration  of  the  frequency  analyzer  was  accomplished  with  a  124  decibel  RAK  pistonphone  prior  to 
testing.  The  pistonphone  produced  a  constant  124  decibel  sound  level  with  a  primary  frequency  of  250  Hz. 
Microphone  calibration  was  checked  periodically  and  no  drift  was  observed.  The  frequency  content  of  thp 
sound  field  was  documented  by  making  filtered  measurements  of  the  sound  level  from  20  to  20,000  Hz  usinq  the 
selective  amplifier  as  a  narrow-band  analyzer.  Frequency  content  of  the  sound  field  was  also  obtained  ny 
sampling  the  microphone  output  directly  and  performing  the  same  analysis  as  described  for  thp  hot-wirp 
signal.  A  comparison  of  the  two  methods  shows  an  excellent  correlation. 

Finally,  to  analyze  the  data  it  was  necessary  to  know  the  frequency  of  the  fan  hladp  passage.  This  was 
obtained  for  each  velocity  and  test  condition  by  marking  a  single  fan  blade  with  a  pipce  of  tapp  and  using  a 
hand-held  strobe  to  determine  the  fan  blade  rpm. 

4. 1.2. 2  Disturbance  Environment  -  no  Airfoil  Present 

Results  of  the  acoustic  and  turbulence  measurements  showed  that  the  experimental  environment  was  a 
complex  function  of  many  variables.  Measurements  were  made  to  take  into  account  as  many  of  these  variables 
as  was  possible^.  Figures  42  and  43  present  the  freestream  turbulence  intensities  and  srund  pressure 
levels  over  the  entire  range  of  velocities.  In  the  velocity  range  from  9  to  31  m/s  no  flow  restrictors  were 
used.  The  turbulence  Intensity  varies  from  0.07%  to  0.151  in  this  range.  Turbulence  intensities  as  low  as 
0.05%  were  calculated  using  analytical  techniques  to  reduce  the  hot-wire  data.  In  using  the  linearlzer  to 
measure  turbulence  intensities  electronic  noise  was  introduced  which  Increased  the  measurements  slightly. 
Therefore  actual  turbulence  intensities  may  be  slightly  lower  than  those  presented  here.  Turhulencp 
intensities  determined  analytically  compared  favorably  with  those  from  the  linearlzer  for  values  greater 
than  0.15%  since  the  electronic  noise  was  no  longer  a  significant  addition  to  the  turbulent  signal. 


There  was  a  noticeable  increase  in  turbulence  intensity  at  a  velocity  of  1?  m/s.  This  corresponds  to  a 
fan  rpm  of  approximately  460.  At  this  fan  rpm  there  was  a  marked  pulsating  of  the  fan  blades  as  the 
belt-drive  from  the  motor  appeared  to  slip  at  this  setting.  The  pulsating  of  the  fan  was  accompanied  hy  a 
slight  squeaking  of  the  belts.  When  a  flow  restrictor  was  introduced  the  turbulence  intensities  in  the 
section  were  increased  significantly.  The  turbulence  intensity  at  idle  speed  (i.e.,  29ft  rpm)  increased  from 
0.07%  to  0.16%  when  one  or  two  flow  restrictors  were  used  (Figure  42).  With  one  flow  restrictor  a  larqe 
increase  in  turbulence  was  observed  at  5.5  m/s  while  an  even  larger  increase  was  observed  at  3. ft  m/s  when 
two  flow  restrictors  were  in  place.  Roth  of  these  velocities  corresponded  to  a  fan  rpm  of  about  460.  It 
may  be  the  increased  work  load  in  conjunction  with  the  pulsating  tendency  of  the  fan  at  this  rpm  which 
causes  the  fan  to  induce  these  high  turbulence  intensities  in  the  test  section.  The  urmb  remains 
essentially  the  same  and  the  decrease  in  U  causes  the  ratio  of  udhc/II  to  increase.  Once  the  460  rpm  region 
is  passed  the  turbulence  intensities  gradually  decreased  hut  always  remain  higher  than  the  no  flow 
restrictor  case.  An  important  comparison  can  be  made  at  0  m/s  (a  chord  Reynolds  number  of  approximately 
150,000  for  the  Lissaman  Airfoil)  where  the  tunnel  may  he  operated  with  no  flow  restrictor  or  with  one  flow 
restrictor  in  place.  The  turbulence  intensity  increases  from  0.07%  to  0.16%  when  one  flow  restrictor  was 
added  after  the  test  section  as  shown  in  Figure  42. 

Analysis  of  the  sound  pressure  levels  in  the  test  section  indicated  a  different  behavior.  Sound  levels 
at  idle  remained  constant  at  approximately  93  dR  (referred  to  2x10*5  N/m?)  regardless  of  the  experimental 
set-up  (Figure  43).  Introduction  of  flow  restrictors  or  a  turbulence  screen  did  not  appear  to  change  the 
total  sound  pressure  level  in  the  test  section  for  a  given  fan  rpm.  However,  to  achieve  the  same  test 
section  velocity  with  one  flow  restrictor  in  place  the  fan  had  to  he  operated  at  a  higher  rpm.  When 
operating  at  a  chord  Reynolds  number  of  150,000  with  one  flow  restrictor  in  place,  a  total  sound  level  of 
104  dB  was  measured  in  the  section  compared  to  a  much  quieter  93  dB  at  the  same  velocity  with  no  flow 
restrictor  in  use.  Special  care  was  taken  to  visually  observe  the  probes  for  mechanical  vibrations  which 
may  have  been  produced  by  trnel  vibrations  or  the  flow  over  the  probe  holders.  The  foam  insulation  between 
the  test  section  and  tunnel  diffuser  successfully  damps  any  mechanical  vibrations  from  the  fan  motor  since 
no  significant  vibration  of  either  the  hot-wire  or  microphone  probes  was  observed. 

Introduction  of  a  turbulence  screen  between  the  test  section  and  tunnel  inlet  increases  the  turbulence 
intensity  in  the  test  section.  Figure  42  shows  the  turbulence  intensities  produced  hy  a  sinqle  turbulence 
screen  with  7.09  meshes/cm  both  with  and  without  one  flow  restrictor  in  place.  The  lowest  set  of  points 
represent  the  case  with  no  turbulence  screens  (i.e.,  the  standard  wind  tunnel  configuration).  At  very  low 
speeds  with  both  the  screen  and  one  flow  restrictor  in  place,  the  screen  did  not  induce  turhulence 
intensities  much  higher  than  those  present  due  to  the  flow  restrictor  alone.  As  the  speed  increased,  the 
turbulence  intensity  increased  significantly.  Turhulence  intensities  produced  hy  the  turhulence  screen  with 
the  flow  restrictor  in  place  were  larger  than  turhulence  intensities  produced  with  the  turhulence  screen 
alone  over  the  same  velocity  range.  This  suggests  that  the  total  turbulence  intensity  may  be  due  to  a 
coupling  of  the  velocity  and  acoustic  fields. 

Total  turbulence  intensity  and  sound  pressure  levels  reveal  many  important  facts  ahout  the  test  'on 
environment.  In  order  to  determine  the  source  of  the  acoustic  and  turbulent  phenomena  in  the  test  se,  •  '.in, 
an  analysis  of  their  frequency  spectra  is  required.  In  addition,  it  is  well  known  that  the  presence  of 
characteristic  frequencies  in  the  freestream  can  affect  the  transition  of  a  laminar  boundary  layer  and 
thereby  radically  change  the  results  of  the  experiment. 

The  frequency  spectra  presented  here  correspond  to  a  chord  Reynolds  number  of  15(1,(10(1  for  the  Lissaman 
Airfoil  (approximately  9  m/s).  Frequency  analysis  was  conducted  with  both  the  sound  pressure  level 
equipment  and  the  hot-wire  anemometer.  Figures  44  and  45  present  frequency  spectra  at  Rc  =  150,000  for  the 
standard  wind  tunnel  configuration  from  the  sound  and  hot-wire  analyses  respectively.  In  these  figures,  as 
well  as  Figures  46  through  47,  the  frequency  spectra  presented  were  normalized  with  respect  to  the  maximum 
value  for  the  case  considered.  Roth  the  sound  and  hot-wire  equipment  picked  up  the  fan  blade  passage 
frequency  of  42  Hz.  This  frequency  was  a  major  part  of  the  acoustic  signal.  While  still  significant  in  the 
turbulent  signal,  it  can  be  seen  that  lower  frequencies  add  a  substantial  amount  of  turhulence  in  addition 
to  the  fan  blade  passage.  These  lower  frequencies  may  include  a  slight  pulsatinq  of  the  fan.  Some 
harmonics  of  the  primary  frequency  can  he  seen  in  both  the  hot-wire  and  acoustic  siqnals.  Figures  46  and  47 
are  frequency  spectra  taken  at  the  same  tunnel  velocity  with  one  flow  restrictor  in  place.  The  fan  blade 
passage  frequency  has  increased  to  ahout  1(10  Hz  as  the  fan  operates  at  a  hiqhpr  spppd  to  compensate  for  the 
pressure  loss  through  the  flow  restrictor.  A  lower  frequency  of  25  Hz  has  also  appeared.  This  frequency  is 
twice  the  fan  rpm  and  becomes  prevalent  due  to  the  slight  pulsating  of  the  fan  under  the  increased  work 
load.  Under  these  test  conditions  the  acoustic  phenomena  have  become  a  larger  part  of  the  turbulent,  signal. 
The  frequencies  associated  with  the  fan  blade  passage  become  the  primary  frequencies  of  the  turbulent  signal 
as  seen  in  Figure  47.  The  lower  frequencies  are  still  present  but  are  of  the  same  or  lower  magnitude  as  the 
fan  blade  frequencies. 

Thus,  part  of  the  increase  in  turbulence  intensity  resulting  from  the  introduction  of  a  flow  restrictor 
is  apparently  due  to  an  increase  in  sound  pressure  level.  Characteristic  frequencies  in  the  freestream  at  a 
given  velocity  (i.e.  Reynolds  number)  will  also  vary  depending  on  whether  a  flow  restrictor  is  used. 

Finally,  with  the  introduction  of  a  turbulence  screen  with  7,09  mpshes/cm  in  the  flow,  the  frequency 
spectrum  becomes  broadband  with  no  characteristic  frequencies  for  both  the  flow  restrictor  and  no  flow 
restrictor  cases  (Figures  48  and  49).  The  turbulence  levels  introduced  by  the  screens  are  much  hiqher  than 
the  free-stream  disturbances  caused  hy  the  fan  blade  passage  or  any  pulsating  of  the  fan.  For  the  Rc  = 
150,000  case  the  turbulence  resulting  from  the  introduction  of  the  screen  dominated  the  test  section 
environment.  This  does  not  mean  the  acoustic  effects  are  absent  but  only  that  the  order  of  this  disturbance 
is  much  lower  than  that  produced  hy  the  turbulence  screen. 
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A  better  understanding  of  the  conditions  present  in  the  test  section  when  flow  restrictors  were 
introduced  to  reduce  the  test  velocity  has  been  obtained.  The  purpose  of  introducing  the  flow  restrictor 
into  the  test  section  was  to  reduce  the  Reynolds  number  of  the  test.  It  was  found  that  some  of  the  earlier 
results  attributed  to  lower  Reynolds  numbers  may  Instead  be  due  to  the  higher  turbulence  intensities  in  the 
test  section.  An  understanding  of  the  small  changes  in  a  test  environment  is  critical  at  these  low  Reynolds 
numbers.  With  the  Introduction  of  flow  restrictors  in  the  test  section  turbulence  intensities  increased 
significantly.  A  part  of  this  increase  could  be  caused  by  the  sound  pressure  waves  transmitted  upstream 
from  the  fan  blade  passage.  The  Increase  in  turbulence  Intensity  at  a  fan  rpm  of  IfiO  was  much  larger  than 
at  any  other  speed.  A  combination  of  Increased  work  load  due  to  the  pressure  drop  through  the  flow 
restrictor  along  with  a  slight  misadjustment  of  the  driv  belts  caused  pulsating  of  the  fan  hlade  and 
subseguent  Increase  in  turbulence.  The  increase  in  turbulence  produced  hy  the  flow  restrictors  had  very 
characteristic  freguencies.  These  frequencies  varied  depending  on  test  conditions.  The  fan  rpm  required  to 
achieve  a  certain  test  section  velocity  was  dependent  on  the  atmospheric  conditions.  The  fan  rpm  also  varied 
with  the  outdoor  wind  velocity  which  impinged  on  the  fan.  Another  factor  which  was  found  to  effect  fan  rpm 
was  small  changes  in  lab  pressure  caused  by  long  per'ods  of  testing  at  high  speeds.  As  the  Notre  name 
tunnels  are  of  the  Indraft,  nonreturn  type,  the  aerospace  laboratory  cannot  he  air  tight.  However,  high 
speed  testing  lowered  the  lab  pressure  enough  to  require  increased  fan  rpm  to  maintain  a  ronstant  velocity. 
This  effect  was  especially  noticeable  with  flow  restrictors  in  place. 

In  contrast  to  the  Increase  in  turbulence  Intensity  produced  hy  the  flow  restrictors,  turbulence 
Increases  due  to  the  introduction  of  a  turbulence  screen  had  no  characteristic  frerjencles  in  the  range 
studied  and  the  intensities  were  relatively  constant  over  the  complete  range  of  velocities.  Turhulence 
produced  by  the  turbulence  screen  dominated  over  that  produced  hy  the  flow  restrictor  when  they  were  used  in 
combination  with  the  latter  except  in  the  case  of  very  low  velocities.  With  the  knowledge  of  these 
experimental  conditions  a  better  understanding  of  data  taken  at  low  Reynolds  numbers  can  he  obtained. 
Further  research  into  the  freestream  turbulence  phenomenon  has  provided  additional  insights  into  the 
character  and  source  of  the  "turbulence"  as  measured  in  the  above  studies.  The  installation  of  improved 
facilities,  most  significantly,  a  DEC  PDP-11/23  computer  based  data  acquisition  system,  has  enabled  the 
refinement  of  measurements  and  data  reduction  schemes  of  similar  experiments  conducted  in  the  University  of 
Notre  Dame  low  speed  wind  tunnels. 

As  stated  earlier,  hot-wire  experiments  have  been  conducted  to  determine  the  level  of  freestream 
turbulence  In  the  testing  environment.  Test  runs  in  which  1024  samples  were  obtained  were  reduced  to 
determine  turbulence  levels  and  frequency  spectra.  It  was  observed  that  a  significant  portion  of  the  energy 
is  tied  up  in  the  lower  frequencies  as  Indicated  by  the  hot-wire  data.  Expansion  of  the  sample  lengths  to 
2048  points  and  the  introduction  of  ensemble  averaging  has  permitted  more  representative  results  to  he 
produced.  Specifically,  there  are  strong  indications  that  for  the  in-draft,  non-return  wind-tunnels,  an 
atmospheric/tunnel  drive  chain  couple  exists  which  tends  to  bias  the  calculated  turhulence  intensities  in  a 
random  fashion.  However,  through  the  use  of  appropriate  filtering  techniques,  this  bias  can  he 
significantly  reduced  and  repeatable  results  obtained.  It  appears  that  for  the  low-speed  Notre  name  wind 
tunnels,  turbulence  intensities  equal  to  or  below  O.lt  are  consistently  obtained  when  a  high-pass  filter  of 
10  hi  is  utilized.  Energy  in  the  frequencies  below  this  value  are  most  likely  due  to  atmospheric 
fluctuations  and  fan  motor  loading  response  to  those  changes  rather  than  to  tunnel  construction  or  design. 
Further,  this  energy  is  several  orders  of  magnitude  larger  than  that  associated  with  frequencies  of  greater 
than  10  Hz.  Thus,  it  is  apparent  that  these  low  frequency  phenomena  do  not  contribute  to  freestream 
turbulence  as  it  is  known  in  the  general  sense  and  so  shoild  not  he  allowed  to  influence  evaluations  of  the 
wind  tunnel  turbulence  level.  However,  it  is  entirely  appropriate  that  when  quoting  the  freestream 
turbulence  level  for  a  specific  wind  tunnel .the  investigators  should  indicate  the  frequency  cutoff  used  in 
its  determination.  Accordingly,  efforts  can  then  concentrate  on  isolating  the  effects  of  the  predominate 
frequencies  in  the  resulting  range  on  more  general  experimental  studies. 

4. 1.2. 3  Airfoil  Performance 

The  lift  and  drag  performance  of  the  smooth  Lissaman  airfoil  in  the  standard  wind  tunnel  configuration 
(i.e.,  no  screen  or  flow  restrictors,  see  Figure  31)  is  shown  in  Figure  50.  As  the  angle  of  attack  was 
Increased,  smoke  visualization  indicated  that  at  fi°  the  laminar  boundary  layer  separated  on  the  upper 
surface  at  about  25%  chord  while  at  8°  the  boundary  layer  appeared  to  he  undergoing  transition  and  separated 
from  the  upper  surface  at  about  35%  chord.  At  an  angle  of  attack  of  10’  transition  appeared  to  he  complete 
and  the  boundary  layer  remained  attached  until  about  the  70%  chord  location.  There  is  a  noticeable  change 
in  the  lift  curve  slope  associated  with  the  extension  of  attached  turbulent  flow.  A  smoke  photograph  at  a  = 
12°  is  shown  in  Figure  51a.  The  lift  coefficient  continues  to  Increase  in  this  region.  Figure  50,  until  it 
reaches  a  maximum  value  of  1.3  at  16°.  Further  increases  in  angle  of  attack  cause  the  location  of  turbulent 
separation  near  the  trailing  edge  to  move  upstream  (and  to  decrease  slightly)  until  it  reaches  about  35% 
chord  where  a  jump  takes  place  to  a  laminar  separation  at  the  leading  edge  at  ahout  14°.  At  this  point 
there  is  an  abrupt  decrease  in  c,  from  1.25  to  ahout  O.b.  As  the  angle  of  attack  is  decreased  from  25°,  the 
boundary  layer  separates  in  the  laminar  state  and  the  ct  remains  ahout  O.R  until  an  angle  of  11°  is  reached. 
With  little  or  no  free  stream  turbulence  present  the  flow  transitions  shortly  after  separation.  A  comparison 
of  the  airfoil  flow  field  it  o-  12°  for  both  Increasing  and  decreasing  angle  of  attack  is  shown  in  Figure 
51.  The  lift  jumps  up  at  a  *  10°  as  a  result  of  the  fact  that  transition  in  the  separated  shear  layer 
allows  the  flow  to  reattach.  The  accompanying  variations  in  the  profile  drag  coefficient  is  shown  in  Figure 
50.  The  abrupt  decrease  in  c,  is  accompanied  by  an  abrupt  Increase  in  c^.  Therefore  in  the  lowest 
turbulence,  quietest  wind  tunnel  configuration,  a  significant  clockwise  or  high  c«nax  hysteresis  region  in 
the  lift  and  drag  forces  was  found.  The  presence  and  extent  of  this  hysteresis  was  determined  by  the 
location  of  separation  and/or  transition  In  the  boundary  layer.  The  location  of  transition  from  laminar  to 
turbulent  flow  in  the  boundary  layer  has  been  known  for  a  long  time^fi  to  be  affected  by  the  level  and  type 
of  free  stream  disturbances. 


In  earlier  experiments  by  Conigliaro^9  using  this  airfoil  ,  hysteresis  was  not  found.  This  data  was 
taken  by  increasing  the  angle  of  attack  from  -10°  to  '•20°  and  then  turning  the  tunnel  off  for  the  balance 
calibration.  The  airfoil  was  then  returned  to  -10°  angle  of  attack  for  the  next  experiment.  In  the  present 
investigation  no  attempt  was  made  to  determine  whether  or  not  hysteresis  occurred  at  negative  angles  of 
attack . 

The  result  of  changing  the  acoustical  environment  by  adding  one  flow  restrictor  at  the  end  of  the  test 
section  is  shown  in  Figure  52.  As  discussed  earlier,  the  addition  of  one  restrictor  increases  both  the  free 
stream  turbulence  level  and  the  sound  pressure  level  for  a  fixed  value  of  tunnel  velocity.  This  test 
section  environment  reduced  the  size  of  the  hysteresis  region  and  produced  a  slightly  higher  c  max  of  almost 
1.4.  A  slightly  lower  minimum  drag  coefficient  was  also  obtained.  The  use  of  two  flow  restrictors  produced 
similar  results  with  the  hysteresis  being  almost  completely  eliminated.  The  increase  in  free  stream 
turbulence  and  acoustic  excitation  caused  the  laminar  shear  layer  to  transition  much  earlier,  thus  allowing 
the  flew  to  reattach  sooner. 

Increasing  the  free  stream  turbulence  level  to  about  0.3%,  by  adding  one  7.09  meshes/cm  screen  at  the 
upstream  end  of  the  test  section  with  no  flow  restrictor,  produced  the  lift  and  drag  coefficients  presented 
in  Figure  53.  This  test  section  environment  completely  eliminated  the  hysteresis  region  and  yielded  values 
of  ctmax  antJ  cdmin  between  those  of  Figures  50  and  52.  With  a  larger  turbulence  intensity  in  the  test 
section,  the  airfoil  boundary  layer  transitions  very  close  to  the  leading  edge,  eliminating  hysteresis  hy 
enabling  the  flow  to  reattach  at  higher  angles  of  attack.  The  abrupt  decrease  in  cooccurred  at 
approximately  the  same  angle  of  attack  in  each  case.  The  very  large  adverse  pressure  gradient  at  this  angle 
of  attack  (i.e.  19°)  caused  the  boundary  layer  to  separate  whether  it  was  laminar  or  turbulent.  Hysteresis 
occurred  because  the  laminar  separated  shear  layer  did  not  reattach.  An  increase  in  turbulence  did  not 
prevent  the  abrupt  loss  of  lift,  but  the  separated  flow  was  turbulent  allowing  more  rapid  reattachment. 

When  the  chord  Reynolds  number  was  increased  to  200,000  the  hysteresis  region  was  reduced  when  using 
the  standard  wind  tunnel  configuration.  At  this  condition  the  abrupt  decrease  in  lift  occurred  at  about  19° 
for  increasing  angle  of  attack  and  the  lift  jumped  up  when  the  angle  of  attack  was  decreased  to  16°.  At  a 
chord  Reynolds  number  of  300,000  the  lift  decreased  abruptly  at  about  21°  and  jumped  hack  up  at  about  ?0°. 

Free-stream  disturbances  are  a  major  source  of  disparity  in  experimental  data.  However,  there  are 
other  sources  of  disparity  which  produce  results  similar  to  those  produced  by  free-stream  turbulence. 
Figures  54  shows  the  lift  and  drag  curves  produced  in  the  standard  wind  tunnel  environment  with  a  strip  of 
tape  2.5  mm  wide  and  0.15  mm  thick  placed  near  the  leading  edge  (i.e  across  the  span  at  1.1%  chord)  of  the 
airfoil.  This  small  boundary  layer  roughness  or  trip  reduced  the  hysteresis  in  a  similar  manner  to  the 
introduction  of  a  flow  restrictor.  The  tape  produces  similar  results  hy  tripping  the  boundary  layer  and 
causing  early  transition.  A  model  with  a  small  amount  of  surface  roughness  nr  Irregularities  in  the  surface 
caused  by  fabrication  defects  could  produce  the  same  results.  The  effect  of  surface  roughness  warrants 
further  study  and  will  be  addressed  later. 

The  section  lift  coefficient  and  sect’on  drag  coefficient  versus  angle  of  attack  for  the  Miley 
airfoil^  for  Rc  =  150,000  and  Rc  =  200,000  are  shown  in  Figums  55  and  56  respectively.  The  lift 
coefficient  Is  zero  at  an  angle  of  attack  of  about  -3°  for  all  of  the  Reynolds  numbers  studied.  At  a  chord 
Reynolds  number  of  150,000  there  is  significant  hysteresis  in  the  lift  and  drag  forces.  The  smoke 
visualization  photographs  shown  in  Figure  57  for  Rc  =  150,000  are  helpful  in  understanding  how  the 
hysteresis  is  produced  in  these  experiments.  At  an  angle  of  attack  nf  zero  degrees,  Figure  57a,  thp 
boundary  layer  separates  near  the  maximum  thickness  location  on  the  upper  surface  while  a  separation  bubble 
near  the  leading  edge  on  the  lower  surface  appears  to  trip  the  flow  on  this  side.  There  is  no  radical 
change  in  the  observed  behavior  of  the  upper  and  lower  surface  boundary  layers  as  the  angle  of  attack 
increases  to  about  9°.  There  is,  however,  some  upstream  movement  of  the  separation  point  and  the  pressure 
distribution  is  altered  to  produce  greater  lift,  etc.  The  hysteresis  region  occurs  between  angles  of  attack 
of  about  10°  and  17°.  Increasing  the  angle  of  attack  to  13°  (Figure  57c),  and  then  to  15°  (Figure  57d) 
moves  the  laminar  separation  point  toward  the  leading  edge  with  very  little  change  in  lift  and  a  continuous 
increase  in  drag.  At  about  17°,  the  lift  increases  dramatically  while  the  drag  decreases.  After  this  jump 
the  lift  drops  off  gradually  while  the  drag  increases.  Decreasing  the  angle  of  attack  from  about  20° 
produces  a  much  larger  lift  and  lower  drag  from  17°  to  10°.  The  reason  for  this  is  evident  in  the  smoke 
visualization  photographs  taken  at  15°  and  13°  when  the  angle  of  attack  was  decreased  (e.g.,  see  Figure 
57 f ) .  For  these  cases  transition  aopears  to  take  place  downstream  of  the  maximum  thickness,  allowing  thp 
upper  surface  boundary  layer  to  remain  attached  almost  to  the  trailing  edge.  Thus  the  airfoil  boundary 
layer  has  an  entirely  different  character  when  the  angle  of  attack  is  increasing  toward  about  ?0°  than  when 
it  is  decreasing  from  about  20°. 

For  chord  Reynolds  numbers  of  200,000  and  higher  (i.e..  Figure  56)  r.o  hysteresis  is  present  and  this 
smooth  airfoil  performs  as  expected.  It  is  clear  from  these  figures  that  operating  the  Miley  airfoil,  which 
was  designed  for  Rc  =  600,000  at  lower  Reynolds  numbers  severely  degrades  Its  performance  and  produces  a 
hysteresis  loop  which  acts  in  the  opposite  sense  as  the  loop  produced  by  the  Lissaman  airfoil. 

Although  the  hysteresis  loop  for  Rc  =  100,000  and  lower  was  documented  and  verified  with  static 
pressure  measurements^4*},  the  surface  roughness  due  to  the  very  small  static  pressure  taps  in  the  pressure 
model  was  enough  to  eliminate  hysteresis  at  Rc  =  150,000.  For  this  reason  the  influence  of  disturbance 
environment  was  also  studied  for  the  Miley  airfoil.  The  addition  of  onp  flow  restrictor  increased  thp  free 
stream  turbulence  level  from  0.07%  to  0.16%  while  the  introduction  of  one  7.09  mesnes/cm  screen  upstream  of 
the  model  raised  the  free  stream  turbulence  from  0.07%  to  about  0.30%.  A  strip  of  tape  0.127  mm  thick  and 
2.21  mm  wide  was  placed  on  the  airfoil  with  the  leading  edge  of  the  tape  at  0.013c.  The  results  indicate 
that  the  hysteresis  loop  is  not  observed.  The  remained  the  •  ame  for  all  configurations  at  Rc  = 
150,000.  It  was  found  than  any  increase  in  the  disturbance  environment  eliminated  thp  hysteresis  loop  at  Rc 
=  150,000. 
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The  importance  of  this  hysteresis  phenomena  cannot  he  overemphasized.  Low  Reynolds  number  airfoil  data 
obtained  in  noisy  and/or  high  turbulence  wind  tunnels  may  not  exhibit  significant  hysteresis.  Therefore, 
aircraft  designed  using  such  wind  tunnel  data  may  not  perform  as  expected  in  fliqht  where  the  free  stream 
disturbance  level  is  often  very  low. 

4. 1.2. 4  Remarks 


The  problems  associated  with  obtaining  accurate  wind  tunnel  data  for  airfoil  sectiuns  at  low  Reynolds 
numbers  are  compounded  by  the  extreme  sensitivity  of  the  boundary  layprs  to  the  free  stream  disturbance 
environment.  The  effect  of  free  stream  disturbances  varies  with  magnitude,  frequency  content,  and  source  of 
the  disturbance.  The  sensitivity  and  accuracy  of  the  measurement  and  data  acquisition  systems  as  well  as 
the  experimental  procedure  used  can  have  a  substantial  effect  on  the  results  obtained. 


Although  free-stream  disturbances  produced  the  largest  disparity  between  different  tests  for  the 
Lissaman  and  Miley  airfoils,  not  all  of  the  differences  can  he  attributed  to  free-stream  disturbances. 
Model  imperfections  or  surface  roughness  can  produce  results  identical  to  those  achieved  with  free  stream 
disturbances.  Reynolds  number  effects  are  critical  at  low  speeds.  An  increase  in  Reynolds  number  from 
150,000  to  200,000  for  the  Lissaman  airfoil  will  eliminate  a  major  portion  of  the  hysteresis,  and  the 
hysteresis  is  insignificant  at  300,000.  For  the  Miley  airfoil  an  increase  in  chord  Rpynolds  number  from 
150,000  to  200,000  completely  eliminates  hysteresis.  It  is  important  that  the  free-stream  disturbances  he 
well  documented  for  each  test  condition  in  order  to  correctly  attribute  differences  in  test  results  to  these 
free  stream  disturbances.  In  accordance  with  new  techniques  and  equipment,  such  disturbances  in  the  Notre 
Dame  tunnels  continue  to  he  investigated  and  recorded.  A  clear  distinction  between  the  effects  nf 
free-stream  disturbances,  model  irregularities,  and  Reynolds  numher  must  he  made  before  the  performance  of 
airfoils  at  low  Reynolds  numbers  can  be  clearly  understood. 

4.1.3  The  Influence  of  Surface  Roughness 


As  mentioned  earlier  the  effects  of  roughness  on  the  performance  of  airfoils  operating  in  the  low 
Reynolds  number  regime  can  be  quite  significant.  However,  to  date  there  has  been  no  systematic  study 
dealing  with  this  problem  and  most  data  that  does  exist  has  only  hern  gathered  as  part  of  general 
evaluations  of  the  performance  of  particular  airfoils*»35>37,43,59t  That  no  emphasis  has  been  placed  on 
investigating  the  effects  of  roughness  in  this  area  is  not  surprising  when  one  considers  that  a  hasic 
understanding  of  the  Influence  of  roughness  on  boundary  layer  development,  stability,  and  transition  is  as 
yet  undeveloped. 


References  141  thru  146  present  a  small  overview  of  some  of  the  more  fundamental  investigations  into 
roughness  effects  on  boundary-layer  behavior.  Unfortunately,  most  work  has  dealt  with  flat  plates  or  pipe 
flow  with  little  attention  given  to  pressure  gradients  over  a  curved  surface.  Tani^l,145  has  done  much  to 
illuminate  the  complexities  of  experimental  studies  into  roughness  and  in  identifying  the  important 
parameters.  Klebanoff  and  Tidstrom*46  have  also  contributed  mucn  in  this  area  including  demonstrating  the 
usefulness  of  applying  stability  theory  to  problems  in  this  field. 


Experimental  data  concerning  the  effect  of  roughness  on  airfoil  performance  has  been  gathered  at  the 
University  of  Notre  Dame  and  elsewhere*4^'*^  Unlike  the  situation  at  highpr  Rpynolds  numbers,  surface 
roughness  can  have  quite  beneficial  effects  on  the  performance  of  airfoils  at  lower  Reynolds  numbers, 
although  this  is  by  no  means  always  the  case.  Typically,  roughness  is  used  to  overcome  the  adverse  pffects 
of  laminar  separation  by  inducing  transition  to  a  turbulent  flow  which  remains  attachpd  to  the  airfoil 
surface.  Pohlen  and  Mueller  observed  elimination  of  the  hysteresis  region  for  the  Miley  MD6-13-12R  airfoil 
at  Rc  =  150,000,  through  the  use  of  a  tape  strip'*0.  In  contrast,  Payne  and  Nelson  recorded  sevpre  losses  in 
the  performance  of  the  NACA  0018  when  at  the  same  Reynolds  number  and  roughened  with  grit*a^.  Carmichael* 
discusses  data  gathered  at  Stuttgart  which  shows  the  effect  of  two-dimensional  tripping  strips  on 
performance  parameters  for  a  variety  of  airfoils.  While  some  of  the  airfoils  studied  showpd  performance 
degradation,  others  demonstrated  remarkable  Improvements. 

Recent  studies  at  the  University  of  Notre  Dame  on  the  Wortmann  FX  63-137  airfoil  have  underscored  snmp 
of  the  complexities  of  roughness  studies  in  the  low  Reynolds  numbpr  regime*^.  The  models  used  had  chord 
lengths  of  152.4  mm  (6  in)  and  spans  of  412.75  mm  (16.25  In).  Tests  conducted  at  a  Reynolds  number  of 
100,000  have  shown  that  roughness  type,  size,  and  location  can  each  have  a  significant  effect  on  airfoil 
behavior.  However,  as  Is  the  case  with  many  other  facets  of  low  Reynolds  numher  testinq,  problems  exist 
with  the  investigative  techniques  themselves  which  must  he  considered.  As  an  example.  Figures  58  and  59  are 
presented  which  show  wide  differences  in  pressure  distribution  measurements  for  two  models  of  the  samp 
airfoil  at  the  same  test  conditions.  One  model  was  constructed  with  pressure  taps  staggered  so  as  to  place 
each  succeeding  one  outside  any  turbulent  wedge  that  might  exist  behind  any  previous  tap.  The  other  model 
has  taps  placed  one  behind  the  other  In  the  chordwlse  direction.  As  seen  in  Figure  59,  the  "in-line"  tap 
model  falls  to  Indicate  the  separation  bubble  which  exists  on  the  FX  63-137  at  16°  angle  of  attack  and  a  Rc 
of  200,000.  Evidentally,  the  taps  of  this  model  act  in  a  similar  fashion  to  surface  roughness  by  causing 
transition  at  a  more  forward  location  than  would  have  been  the  case  had  a  bubble  formed  and  tripped  the  flow 
further  aft. 


Additional  pressure  plots  shown  In  Figures  60  and  61  demonstrate  the  effect  of  surface  grit  on  the 
pressure  distribution  of  the  FX  63-137  at  16°  and  Rc  =  100,000  (k  is  the  height  of  the  grit  and  w  is  the 
width  of  the  grit  band).  As  seen  here,  grit  of  sufficient  height  can  cause  stall  at  this  angle  in  contrast 
to  the  clean  airfoil  which  exhibits  stall  at  19°.  Furthermore,  even  for  the  smallest  qrit  size  shown, 
transition  is  moved  forward  from  the  0.18c  location  seen  on  the  smooth  airfoil  to  approximately  0.09c  on  the 
roughened  model. 
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In  the  course  of  these  studies,  the  Notre  Dame  force  balance  was  improved  by  the  addition  of  a 
moment-measuring  capability  which  utilized  strain  gages.  Figures  62  and  63  provide  force  halance  data  for 
the  FX  63-137  in  the  clean  and  rncghened  configurations.  In  the  latter  case,  a  strip  of  tape  composed  of 
several  layers  of  tape  was  notched  to  produce  a  zig-zag  pattern  and  was  placed  across  the  upper  airfoil 
surface  at  x/c  =  0.01.  In  this  comparison,  C?max  is  significantly  decreased  and  (I^min  slighvly  increased. 
However  ( £/d ) ma x  is  increased  by  over  7%  although  it  occurs  at  6°  instead  of  the  8°  observed  for  the  case  of 
the  smooth  airfoil.  Differences  in  the  variation  of  quarter  chord  moment  coefficient  with  angle  of  attack 
can  also  be  seen.  The  angle  of  attack  region  of  hysteresis  is  reduced  to  just  two  degrees  whereas  for  the 
smooth  airfoil  it  occurred  across  four  degrees. 

Given  the  complexity  of  the  low  Reynolds  number  problem,  it  does  not  appear  tnat  general  solutions  for 
the  prediction  of  airfoil  performance  under  roughened  conditions  will  be  available  for  some  time.  Instead, 
a  more  basic  understanding  of  the  influence  of  roughness  on  boundary-layer  behavior  will  have  to  be 
developed.  Until  such  an  understanding  is  obtained,  the  design  engineer  working  in  the  field  of  low 
Reynolds  number  airfoils  will  have  to  depend  on  wind  tunnel  results  of  the  type  discussed  above. 

4.1.4  Finite  Wing  Experiments 

In  order  to  design  a  practical  low  Reynolds  numher  vehicle,  It  is  necessary  to  evaluate  the  importance 
of  three-dimensional  boundary  iayer  effects  on  a  finite  wing.  Many  of  t..ese  vehicles  require  larqp  aspect, 
ratio  wings;  however,  the  largest  aspect  ratio  allowable  in  most  wind  tunnels,  usinq  a  reasonable  chord 
jimension,  is  usually  much  smaller  than  the  actual  vehicle.  These  small  aspect  ratio  wing  experiments  can, 
nevertheless,  illuminate  the  major  problems  with  the  finite  wing. 

Studies  have  been  in  progress  at  the  University  of  Notre  Dame  related  to  the  performance  of  finite, 
rectangular  wings^5i  involving  the  Wortmann  FX  63-137  airfoil  section  at  chord  Reynolds  numhers  less  than 
200,000.  Attempts  have  been  made  to  determine  the  effect  of  varying  winq  aspect  ratio  (based  on  semi-spanl 
and  chord  Reynolds  number  on  the  performance  of  the  wing  as  well  as  the  effect  on  t.he  separation  bubble  that 
forms  under  these  conditions.  It  was  first  necessary  to  determine  the  two-dimensional  characteristics  of 
the  profile  to  make  adequate  assessments  of  finite  wing  performance. 

Force  balance  data  was  obtained  for  chord  Reynolds  numbers  between  110,000  and  200,000  for  three  winqs 
with  different  aspect  ratios  ,  as  well  as  tor  the  two-dimensional  case.  No  wind  tunnel  corrections  were 
applied  to  the  force  data  for  the  finite  wing  cases.  Figures  64-66  show  the  two-dimensional  lift,  draq  and 
quarter-chord  moment  coefficients  versus  angle  of  attack  as  well  as  the  drag  polar  for  chord  Reynolds 
numbers  of  80,000,  100,000  and  200,000.  Figure  67  presents  some  of  these  characteristics  at  Rc  =  100,000 
for  a  finite  wing  with  semi-span  aspect  ratio  of  2.7.  From  data  of  this  type,  it  was  possihle  to  identify 
the  effect  of  Reynolds  number  and  aspect  ratio  upon  various  performance  parameters. 

Figure  68  shows  the  lift  and  drag  curves  for  the  infinite  ( i . e . ,  two-dimensional  easel  and  three  finite 
wing  cases  at  Rc  =  200,000.  It  is  evident  that,  as  aspect  ratio  decreases,  the  slope  of  the  lift  curve 
decreases,  the  maximum  lift  decreases,  minimum  drag  increases,  while  the  stall  anqle  of  attack  increases  and 
the  range  of  the  hysteresis  loop  increases. 

It  was  noticed  that  the  performance  of  the  lower  aspect  ratio  winqs  was  a  function  of  chord  Reynolds 
number.  Figure  69  shows  the  lift  and  drag  curves  for  a  finite  winq  of  spmi-span  aspect  ratio  of  two,  at 
chord  Reynolds  numbers  of  80,000,  100,000,  160,000  and  200,000.  The  effect  of  increasing  Reynolds  number 
for  a  given  wing  was  similar  to  the  effect  of  increasing  asppet  ratio  for  a  given  Reynolds  number;  maximum 
lift  increases,  minimum  drag  decreases  while  the  angle  of  attack  of  stall  increases  and  the  range  of  the 
hysteresis  loop  increases. 

The  effect  of  aspect  ratio  and  Reynolds  number  upon  the  maximum  lift  coefficient  is  shown  in  Figurp  70. 
Also  plotted  is  the  two-dimensional  data  from  Reference!-^,  which  compares  reasonably  well.  Increases  in 
Reynolds  number  and  aspect  ratio  causes  a  rise  in  dLmax>  particularly  at  the  lower  aspect  ratios  and 
Reynolds  numbers.  Figure  71  is  a  similar  plot  for  the  minimum  drag  coefficient.  As  was  expected,  increases 
in  Reynolds  number  and  aspect  ratio  decrease  this  parameter.  A  much  lower  minimum  drag  coefficient  is  found 
in  the  Stuttgart  experiments  using  the  wake  rake  technique^-*.  This  data  is  also  shown  in  Fiqure  71.  In 
both  cases  the  greatest  loss  in  performance  came  between  the  largest  aspect  ratio  finite  wing  and  the 
infinite  wing  case. 

Preliminary  data  indicted  that  the  quarter-chord  moment  coefficient  was  not  drastically  affected  by 
the  downwash  field  general,  u  on  finite  wings.  Comparison  of  Figures  65  and  67  show  that  the  quarter-chord 
moment  coefficient  versus  angle  of  attack  curve  became  "smoother"  for  the  finite  winq,  hut  was  not  alterpd 
significantly  in  magnitude. 

Further  studies  of  this  problem  indicated  that  performance  of  the  airfoil  and  the  winqs  was  linked  to 
the  development  of  a  separation  bubble  on  the  upper  surface.  Since  force  halance  data  yielded  very  little 
information  on  this  phenomenon,  other  experiments  were  necessary.  Using  a  specially  desiqned  and 
constructed  pressure  tap  model,  described  in  detail  in  Appendix  A  of  Reference  166,  chordwise  pressure 
distributions  were  obtained  at  seven  spanwise  stations  for  a  finite  wing  (AR=2.D),  as  well  as  for  the 
two-dimensional  airfoil,  at  twojleynolds  numhers.  Figures  7?  and  73  show  the  pressure  distributions  alnnq 
the  wing  span  on  a  finite  wing  {’AR=2 .0) ,  as  well  as  the  two-dimensional  distributions,  at  Reynolds  numbers 
of  80,000  and  200,000.  It  should  be  pointed  out  that  the  static  pressure  taps  of  this  model  were  aliqned 
parallel  to  the  chordline,  hence  an  effective  roughness  was  introduced  into  the  problem.  However,  since  all 
pressure  distributions  in  Figures  72  and  73  were  taken  with  the  same  model,  comparisons  of  this  data  snujld 
yield  valid  information  concerning  separation  bubble  trends.  Laminar  separation,  free  shear  layer 
transition,  turbulent  reattachment  and  turbulent  separation,  if  they  existed,  are  marked  on  the  pressure 
distributions  with  "S",  "T",  "R"  and  "TS"  respectively. 


Figure  72  shows  the  effect  of  the  downwash  field  upon  the  sectional  wing  characteristics  for  Rc  = 
80,000  and  an  angle  of  attack  of  18°.  One  would  expect  the  Induction  of  air  by  the  tip  vortex  to  reduce  the 
effect'  ?  angle  of  attack,  with  the  greatest  reduction  occurring  near  the  tip.  Associated  with  this 
occurrence,  one  would  expect  a  reduction  In  the  suction  peak  as  well  and  laminar  separation,  transition  and 
reattachme.it  occurring  farther  downstream  on  the  model,  as  one  progresses  along  the  span,  from  the  root  to 
the  tip.  The  lower  left  plot  Is  the  farthest  outboard  station,  the  span  station  moving  Inhoard  toward  the 
root  as  one  progresses  up  and  to  the  right  on  the  figure,  until  reaching  the  two-dimensional  distributions 
In  the  upper  right.  The  effect  on  the  pressure  distribution  Is  small  from  root  to  midspan.  All 
characteristics  expected  are  visible  In  Figure  72.  A  similar  plot  Is  presented  In  Figure  73,  for  Rc  » 
200,000  and  an  angle  of  attack  of  18°.  Again,  reduction  of  the  effective  angle  of  attack  was  apparent  In 
the  pressure  distributions  as  the  tip  was  approached.  In  this  case  however,  the  Indications  of  the 
separation  bubble  In  the  pressure  distributions  were  masked  due  to  the  pressure  tap  locations  and 
pecularltles  of  bubble  formation  at  this  condition.  These  distributions  were  integrated  n  yield  the 
sectional  lift  coefficients  at  those  stations.  A  summary  of  all  angles  of  attack  tested  for  Rc  »  2fin,nnn  is 
presented  In  Figure  74.  One  can  see  very  little  load  variation  until  the  fin*  spanwise  position. 

In  order  to  verify  the  flow  pattern  variations  along  the  span  of  the  finite  wings,  surface  oil 
visualization  was  performed.  Lcvlng  and  Katzoff*66,  recommended  that  kerosene  he  used  for  flow  speeds  less 
than  35  m/s,  but  preliminary  studies  indicated  that  this  method  was  not  practical  for  the  low  flow  speeds 
employed  In  this  study  (9  and  19  m/s).  W1nkelmann*6'•*6',  indicated  that  visualization  with  other  thin  oils 
was  possible  at  low  speeds.  A  mixture  of  propylene  glycol  and  water  was  used.  A  fluorescent  dye,  Hranine, 
was  added  as  a  marker,  as  was  a  small  guantlty  of  a  wetting  agent  to  reduce  surface  tension  In  the  fluid. 
This  method  worked  reasonably  well,  except  at  the  higher  angles  of  attack  and  lower  speeds.  Spanwise 
variations  of  flow  patterns  were  gulte  evident  with  this  method  and  general  trends  were  identified.  Data 
obtained  with  this  method  yielded  similar  separation  bubble  fluid  patterns  to  that  of  References  157  and 
158,  although  at  much  lower  speeds.  It  was  concluded  that  this  method  could  be  a  valuable  tool,  provided 
that  it  was  used  with  care. 

4.2  Flight  Research  Experiments 

Because  of  the  extreme  sensitivity  of  low  Reynolds  number  airfoil  boundary  layers  to  the  disturbance 
environment  in  wind  tunnels,  evaluation  of  such  airfoils  in  the  atmosphere  Is  appealing.  Attaching  the  test 
airfoil  section  and  the  necessary  Instrumentation  to  a  land  vehicle  (l.e.,  automobile  or  truck)  or  an 
airborne  vehicle  (l.e.,  sailplane  or  special  RPV)  are  the  obvious  choices.  Land  vehicle  testing  raises  the 
question  of  the  possible  vibrational  Inputs  to  the  model  being  tested.  Sailplanes,  on  the  other  hanrt.  have 
been  used  successfully  for  many  years  as  a  test  bed  for  low  speed  aerodynamics  research*36.*6".*60.*6* . 
Powered  (l.e.,  human2*  or  internal  combustion  engine6)  aircraft  have  been  used  to  determine  the  aerodynamic 
as  well  as  stability  and  control  characteristics  of  low  speed  vechlcles.  Recently,  a  small  light-weight 
powered  RPV  has  been  proposed  as  a  test-bed  for  obtaining  low  Reynolds  number  airfoil  data*'. 

In  order  to  evaluate  a  new  airfoil  design  rather  than  an  entire  vehicle.  It  Is  necessary  to  Isolate  the 
airfoil  model  from  the  test-bed  vehicle  and  to  "fly"  the  airfoil  model  at  the  desired  Incidence  and  chord 
Reynolds  number.  Successful  flight  tests  of  this  type  were  performed  by  Miiey*36  to  evaluate  his  «Ofi-l3-l?R 
airfoil  design.  The  vehicle  used  was  a  Schwelzer  TG-3A  sailplane.  This  glider  has  an  operating  equivalent 
airspeed  range  from  72.4  to  144.8  km/h  (45  to  90  mph)  and  could  be  controlled  to  within  0.40  km/hr  (0.25 
mph)  by  an  experienced  pilot.  After  considering  the  best  location  of  the  test  airfoil  together  with  the 
effect  of  this  configuration  on  the  flight  characteristics  of  the  glider,  Mlley  located  the  40fi  mm  chord 
1270  mn  span  airfoil  model  upstream  and  above  one  of  the  sailplane's  main  wings  at  a  sufficient  distance 
from  the  tutelage  to  eliminate  Interference  effects.  The  test  airfoil  had  fi09.fi  mm  square  side  plates  on 
each  end  to  help  simulate  the  Infinite  airfoil  case.  Static  pressure  data  were  obtained  from  the  40  static 
pressure  taps  In  the  model  and  transition  locations  were  obtained  using  a  moveable  hot-wire  anpmometer 
probe.  All  of  the  necessary  Instrumentation  was  carried  on-board  the  sailplane. 

Flight  data  were  obtained  for  a  chord  Reynolds  number  of  600,000  (the  design  value  for  the  M0fi-13-12R 
airfoil)  at  six  different  angles  of  attack.  A  comparison  of  one  of  the  pressure  distributions  obtained  in 
these  sailplane  tests  and  data  obtained  In  a  wind  tunnel  at  the  University  of  Notre  name  with  the 
theoretical  values  Is  shown  In  Figure  751,10.  Recause  this  case  represents  a  speed  slightly  beyond  the 
normal  upper  limit  of  the  Notre  Oame  low  speed  wind  tunnel,  the  wind  tunnel  pressure  coefficients  are 
Included  only  to  indicate  the  slmi’ar  trend  with  the  flight  data. 

To  obtain  airfoil  data  on  reasonably  sized  models  (l.e.,  a  minimum  of  about  152  mm  chord)  and  Reynolds 
numbers  below  500,000,  piloted  sailplanes  with  good  low  speed  stability  and  control  characteristics  must  he 
used.  Also,  remotely  piloted,  considerably  smaller  sailplanes  may  be  used.  As  the  size  of  the  sailplane 
decreases,  the  ability  to  carry  pressure  measuring  Instrumentation,  etc.  for  airfoil  measurements  Is 
decreased.  For  small  sailplane  or  powered  RPV  test-bed  configurations,  the  desired  airfoil  characteristics 
must  then  be  backed-out  of  the  entire  vehicle  characteristics  from  (for  example)  power-off  sink  rate 
experiments.  Although  this  has  been  accomplished  with  some  success  at  chord  Reynolds  numhers  ahove  fintl.non, 
further  studies  are  necessary  to  demonstrate  the  validity  of  this  approach  at  lower  Reynolds  numhers.  The 
question  of  accuracy  of  data  obtained  and  possible  Interference  from  other  components  of  the  vehicle  must, 
also  be  considered. 

5.  PRACTICAL  APPLICATIONS 

A  wide  variety  of  practical  applications,  where  low  Reynolds  number  effects  are  present  for  all  or  part 
of  the  operating  period,  have  been  studied.  The  success  of  many  of  these  configurations  is  stronqly 
dependent  upon  achieving  the  highest  aerodynamic  efficiency  of  the  lifting  surfaces  l.e.,  the  wings  or 
rotors.  A  few  examples  of  specific  types  of  configurations  and  their  missions  will  be  presented  to 
Illustrate  this  point. 


5.1  Remotely  Piloted  Vehicles  and  Sailplanes 


An  extremely  large  number  of  RPV's  have  been  designed  and  flown  over  the  past  55  years.  Many  of  these 
vehicles  are  little  more  than  small  scale  airplanes  and  are  often  used  to  study  the  flight  performance  and 
stability  characteristics  of  the  full  size  airplane.  Tethered  RPV's  in  large  wind  tunnels  and  untethered 
RPV's  in  vertical  spin  tunnels  and  in  atmospheric  flight  have  been  operated  successfully  for  many  years. 
Recent  advances  in  reducing  the  size  and  weight  of  electronic  payloads  while  increasing  their  sensing 
capability,  the  use  of  light-weight  and  composite  structures,  and  in  development  of  non-conventional 
propulsion  systems  (e.g.,  solar/electric,  fuel  cell,  and  microwave/electric)  have  greatly  enlarged  the 
possible  applications  of  RPV's. 

Preliminary  design  studies  of  high  altitude  aircraft  platforms  (HAAPs)  have  recently  been 
completed111*  14,162-164^  One  possible  mission  for  such  a  solar  powered  vehicle  is  agricultural  surveilance. 
For  this  type  of  mission  a  high  altitude  long  endurance  platform  capable  of  carrying  a  112.5  kg  (250  lbm) 
payload  with  a  daytime  power  requirement  of  300  watts  has  been  studied15?  an(j  an  eariy  design  resulting  from 
this  continuing  effort  is  shown  in  Figure  76.  This  platform  would  be  capable  of  maintaining  20  km  (65,600 
ft)  altitude  over,  for  example,  California's  San  Joaquin  Valley  (32°  -3fiuN  latitude)  for  up  to  12  months. 
During  this  period  the  payload  would  monitor  crop  conditions  on  an  hourly  basis.  Design  studies  indicate 
that  such  a  vehicle  might  have  a  wing  span  on  the  order  of  100  m  ( 328  ft)  and  a  mean  aerodynamic  chord  on 
the  order  of  3  m  (9.84  ft).  At  speeds  of  10  m/s  (32.8  fps)  and  at  20  km,  (65,600  ft)  the  chord  Reynolds 
numbers  are  less  than  500,000.  For  lack  of  suitable  low  Reynolds  number  airfoil  designs,  these  studies  have 
used  published  airfoils  designed  for  other  purposes,  e.g.  the  Liebeck  L  1003M15?  t  the  wortmann 
FX74-CL5-4015?  and  the  Wortmann  FX  63-137164.  Not  only  have  these  airfoils  not  been  designed  for  this  type 
of  vehicle,  but  sufficient  wind  tunnel  and  flight  test  data  are  not  available  to  guarantee  the  actual 
performance  of  these  designs.  Even  with  future  Improvements  in  regenerative  solar  propulsion  systems,  it 
appears  that  more  efficient  airfoil  and  wing  designs  are  necessary.  HAAP  lift  to  drag  ratios  on  the  order 
of  50-60  appear  to  be  necessary.  The  highest  aerodynamic  efficiency  of  any  man-made  vehicle  Is  the  modern 
sailplane  where  lift  to  drag  ratios  of  near  50  have  already  been  attained.  In  order  to  obtain  further 
improvements  in  L/D,  it  appears  that  active  boundary  layer  control  (e.g.,  blowing  or  suction)  and  the  use  of 
multi-element  airfoil  geometries  may  be  necessary.  Furthermore,  the  propeller  used  in  such  HAAPs  have  a 
mean  aerodynamic  chord  which  Is  at  least  an  order  of  magnitude  smaller  than  the  wing  (i.e.,  ahout  0.3m). 
Rotating  at  between  100  and  200  rpm,  most  of  the  propeller  designs  operate  continuously  in  the  low  Reynolds 
number  regime.  Although  there  has  been  some  recent  attention  to  the  design  of  propellers  of  this  type1  , 
the  airfoil  sections  are  also  chosen  from  published  data  and  not  designed  specifically  for  the  new 
application.  The  design  of  this  type  of  propeller  is  similar  to  the  design  of  low  solidity  wind  turhine 
blades  and  will  be  discussed  in  the  next  section. 

The  Mini-Sniffer5  shown  in  Figure  77  is  an  example  of  an  operationally  successful  mini-RPV.  The 
original  mission  called  for  a  vehicle  to  carry  a  26.8  kg  (59  lbm)  air  sampler  to  21,336  m  (7D,nno  ft)  and 
cruise  there  for  an  hour  within  a  range  of  311.8  km  (200  miles)  or  glide  hack  after  a  27.4  km  (PD.OOD  ft) 
climb.  Using  a  hydrazine  monopropellant  engine,  a  76  kg  (167  lbm)  vehicle  was  developed  using  the  Miley 
M06-1 3-128  airfoil  section155  for  the  main  wing  [mean  aerodynamic  chord  =  0.661  m  (2.17  ft),  span  =  5.4R  m 
(18  ft)]  and  the  Clark-V  airfoil  section  for  the  2  bladed  propeller.  The  chord  Reynolds  numbers  for  such  a 
mission  could  get  as  low  as  about  60,000.  In  the  test  program,  however,  [in  which  altitudes  of  about  6. OR 
km  (20,000  ft)  were  attained]  the  minimum  chord  Reynolds  number  was  about  250,000.  The  maximum  L/D  of  this 
vehicle  was  determined  from  power-off  sink  rate  tests  from  6.09  km  (20,000  ft)  to  be  about  15.  (L/D)max 
with  the  landing  gear  off  was  estimated  to  be  slightly  over  20.  Current  mini-RPV  requirements  indicate  that 
improvements  in  aerodynamic  efficiency  in  the  vicinity  of  C5/?/C[)  »  20  are  necessary  to  Increase  on-station 
endurance  of  reciprocating  engine/propeller  vehicles  to  ahout  10  hours.  Furthermore  these  vehicles  must  be 
capable  of  flying  in  all  kinds  of  weather  without  significant  reductions  in  performance. 

5.2  Wind  Turbines 

Over  the  past  ten  years,  development  of  wind  energy  systems  has  been  increasing.  This  growth  is,  in 
part,  due  the  the  Federal  Wind  Energy  Program  in  the  United  States  and  an  interest  in  renewable  energy 
sources  caused  by  the  rising  cost  of  fossil  fuels  worldwide.  The  focus  of  this  work  has  been  on  high 
efficiency,  lift-type  wind  turbines  such  as  the  HAWT  (horizontal  axis  wind  turhine)  and  the  VAWT  (vertical 
axis  wind  turbine)  shown  in  Figure  78.  Roth  types  are  designed  to  extract  energy  through  an  efficient 
aerodynamic  interaction  between  the  wind  and  blade  elements.  This  optimization  has  led  to  a  class  of 
turbines  characterized  by  high  efficiency,  high  speed  (i.e.,  tip  to  wind  speed  ratios  of  ahout  6,  shown  in 
Figure  79),  and  low  solidity  (ratio  of  blade  area  to  swept  area).  These  characteristics  of  the  modern  wind 
turbine  lead  to  several  consequences. 

Until  now,  the  structural  integrity  of  these  wind  machines  has  been  the  primary  focus  of  research  due 
to  the  high  blade  loadings.  Here  the  aerodynamics  were  considered  only  as  input  to  aeroelastic  problems. 
Choice  of  airfoil  section  was  generally  based  upon  structural  considerations  and  the  availability  of  airfoil 
section  data.  This  data  was  used  not  only  in  structural  analysis,  but  also  in  predicting  turbine 
performance.  Unfortunately,  most  available  section  data  was  for  Reynolds  numbers  of  several  million.  In 
both  the  HAWT  and  VAWT,  chord  Reynolds  numbers  of  less  than  500,000  do  occur  in  sections  of  the  rotor  near 
the  rotating  axis.  Modern  turbines  with  power  output  less  than  about  15  kW  may  experience  Reynolds  numbers 
less  than  500,000  at  all  times  due  to  their  low  solidity.  The  choice  of  a  low  Reynolds  number  airfoil  for  a 
wind  turbine  involves  several  considerations  which  are  unique  to  this  application. 

The  problems  associated  with  the  design  of  a  low  Reynolds  number  HAWT  are  the  same  as  those  encountered 
in  designing  a  low  Reynolds  number  propeller.  Essentially,  the  design  of  the  airfoil  section  must  take  into 
account  the  operating  character  of  the  rotor.  Lift  and  drag  produced  by  a  HAWT  rotor  blade  element  may  be 
resolved  into  an  equivalent  normal  and  thrust  force  as  shown  in  Figure  80.  The  magnitude  of  the  driving 
force  is  a  function  of  the  relative  velocity  and  angle  of  attack.  Additionally,  one  must  consider  the 
geometric  twist  and  relative  velocity  as  'unctions  of  the  blade  element  radius.  Airfoil  design  is 


complicated  in  two  ways.  First,  the  Reynolds  number  is  a  function  of  radius.  This  is  equivalent  to  a 
finite  wing  experiencing  spanwise  wind  shear.  Second,  the  direction  and  speed  of  the  wind  are  not  constant, 
causing  an  unsteady  loading  of  the  turbine.  Also,  consideration  must  be  given  to  three  dimensional  effects 
due  to  the  finite  aspect  ratio  and  the  geometric  twist  of  the  blade. 

The  Darrieus  VAWT  has  more  complex  aerodynamics.  This  type  of  wind  turbine  produces  power  as  a  result 
of  the  tangential  thrust  as  shown  in  Figure  81.  Operation  is  complicated  by  the  fluctuating  angle  of  attack 
experienced  by  a  blade  element.  The  designer  must  not  only  account  for  the  design  state  of  the  wind 
turbine,  but  also  an  unsteady  pitching.  In  order  for  power  to  be  produced  at  all  azimuthal  blade  angles, 
the  thrust  force  must  be  positive  for  both  positive  and  negative  angles  of  attack,  hence,  one  would  expect  a 
symmetrical  airfoil  to  be  practical.  Again,  for  structural  considerations  and  data  availability,  NAOA 
symmetrical  sections  are  currently  used.  Some  testing  has  been  done  on  these  airfoils  at  low  Reynolds 
numbers*®®,  but  little  has  been  done  to  address  the  unsteady  behavior  of  these  airfoils  at  low  Reynolds 
numbers.  In  addition  to  the  unsteady  angle  of  attack  phenomenon,  the  range  of  angle  of  attack  may  approach 
180°  in  portions  of  the  rotor  near  the  axis.  Section  data  for  full-range  (-18(1°  *  18(1°)  is  very  limited 
even  for  high  Reynolds  numbers.  Example  of  thrust  coefficient  data  versus  angle  of  attack  for  a  symmetrical 
NACA  0015  airfoil  at  Rc  =  350, 000*®?  and  for  a  highly  cambered  Wortmann  FX  63-137  airfoil  at  Rc  =  125,000168 
are  shown  in  Figure  82.  The  FX  63-137  data  was  obtained  by  changing  the  anqle  of  attack  only  in  one 
direction,  therefore,  no  hysteresis  could  be  observed. 

The  design  and  testing  of  low  Reynolds  numher  sections  for  use  in  wind  turbines  requires  careful  study 
of  the  rotor  aerodynamic  characteristics.  In  general,  the  desiqner  must  consider  unsteady, 

three-dimensional  performance  over  a  wide  range  of  Reynolds  numbers  and  angles  of  attack.  There  is  a 
definite  need  for  full-range  low  Reynolds  numher  airfoil  section  data  for  use  in  wind  turbine  desiqn  codes. 
Data  is  also  needed  for  full-range  unsteady  airfoil  section  characteristics.  To  this  end,  a  study  is 
currently  in  progress  at  the  University  of  Notre  Dame  addressing  the  characteristics  of  airfoils  at  low 
Reynolds  numbers  in  unsteady  flows. 

6.  CONCLUDING  REMARKS 

In  order  to  achieve  the  aerodynamic  efficiency  levels  necessary  for  successful  operation  o*  low 
Reynolds  number  vehicles  at  subcritical  speeds,  improved  airfoil  sections  and  wing  desiqns  are  necessary. 
Present  design  techniques  produce  airfoil  sections  which  perform  as  they  are  designed  for  chord  Reynolds 
numbers  above  about  500,000.  For  chord  Reynolds  numbers  be’ow  anout  SOO.ODO,  airfoil  sections  are  much  morp 
difficult  to  design  with  confidence  because  the  effects  of  laminar  separation  huhhles  and/or  transition  have 
not  been  adequately  modelled  in  the  design  and  analysis  methods.  The  fact  that  more  than  one  airfoil 
flowfield  is  possible  at  a  given  angle  of  attack  {f.e.,  hysteresisl  indicates  that  further  research  is 
necessary  in  this  area. 

To  improve  the  design  and  analysis  methods  in  this  flight  regime,  experimental  determination  of  thp 
influence  of  free  disturbance  environment,  surface  roughness,  and  unsteady  flow  effects  is  necessary.  Dur 
understanding  of  these  complex  flow  phenomena  can  he  enhanced  with  definitive  wind  tunnel  experiments.  In 
wind  tunnel  experiments,  care  must  be  exercised  to  determine  the  influence  of  the  test  section  disturbance 
environment,  the  experimental  technique,  instrumentation,  data  processing  and  procedure  used,  on  the 
results.  The  airfoil  boundary  layer  may  also  be  sensitive  to  small  imperfections  in  the  wind  tunnel  model 
profile.  Flight  experiments,  though  usually  more  time  consuminq  and  expensive  than  wind  tunnel  experiments, 
may  also  be  necessary  to  insure  that  a  given  airfoil  dpsign  receives  a  fair  evaluation  under  actual 
operating  conditions. 

7.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

Given  the  wide  variety  of  low  Reynolds  number  applications  of  current  interest  and  thp  present  small 
aerodynamic  data  base,  the  following  recommendations  are  presented  for  future  research.  Roth  single-arid 
multi-element  airfoil  sections  need  to  be  considered  as  well  as  the  benefits  nf  in  Ive  boundary  layer 
control  concepts. 

1)  The  question  of  when  and  where  transition  from  laminar  to 
turbulent  flow  takes  place  in  the  attached  or  separated 
airfoil  boundary  layer  Is  of  primary  importance. 

Improvements  in  both  design  and  analysis  methods  require  a 
more  complete  understanding  of  the  stability  and  transition 
of  viscous  shear  layers.  Definitive  experiments  into  these 
phenomena  in  airfoil  type  pressure  gradients  as  affected  by 
the  free  stream  disturbance  environment  and  airfoil  surface 
roughness  characteristics  are  necessary.  All  wind  tunnel 
experiments  in  this  area  must  be  preceeried  by  an  extensive 
study  of  the  unique  disturbance  environment  of  the 
experimental  facility  used  and  how  this  may  affect  the 
results. 

2)  In  order  to  include  a  separation  bubble  model  in  both  design 
and  analysis  methods,  more  detailed  experimental  studies  are 
needed  of  the  structure  of  the  rather  extensive  separation 
bubbles  which  occur. 

3)  A  systematic  study,  both  experimental  and  analytical,  of  the 
unsteady  flow  effects  on  laminar  separation  and  transition 
is  also  necessary. 


4)  To  improve  the  experimental  evaluation  of  new  airfoil  and 
vehicle  designs,  the  experimental  measurement  techniques 
should  continue  to  be  improved  and  evaluated  so  that  more 
accurate  wind  tunnel  and  flight  data  may  be  obtained. 

5)  Analytical  methods  should  continue  to  be  extended  and 
improved  in  conjunction  with  experimental  findings.  The 
analytical  and  experimental  approaches  should  proceed  so  as 
to  compliment  each  other.  In  this  way  a  more  rapid  solution 
to  the  design  problems  will  be  achieved. 

6)  The  analytical  and  experimental  techniques  should  be 
extended  to  include  the  fully  integrated  practical 
configurations  of  interest. 
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Figure  1.  Chord  Reynolds  Number  versus  Flight  Velocity  for  a  Variety  of 

Natural  and  Man-Made  Flying  Objects,  Adapted  From  References  3  and 
4.  (All  values  based  on  standard  sea  level  conditions  except  the 
high  altitude  aircraft  platforms  where  values  were  based  on 
conditions  at  about  70,000  feet). 
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Figure  2.  Maximum  Aerodynamic  Efficiency  versus  Speed  for  a  Variety  of 

Natural  and  Man-Made  Flying  Objects  -  Adapted  from  Reference  4. 
(The  Shaded  Regions  Represent  Estimates  for  HAAPs  and  Minl-RPVs). 
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Figure  3.  Wing  Area  Versus  Mass  for  a  Variety  of  Natural  and  Man-Made  Flying 
Objects  (Reference  30). 
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Figure  4.  High  Reynolds  Number  Flow  Over  a 
Large  Aircraft  Wing. 
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Figure  5.  Low  Reynolds  Number  Flow  Over  a 

Small  RPV  Wing  at  a  Small  Angle  of 
Attack . 
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Figure  6.  Low  Reynolds  Number  Flow  Over  a 

Small  RPV  Wing  at  a  Large  Angle  of 
Attack  (Chordwise  Length  of 
Separation  Bubble  Greatly 
Exaggerated) . 


Figure  7.  Low  Reynolds  Number  Flow  Over 
Small  RPV  Wing  at  Stal  1 . 
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Figure  8.  Effect  of  Long  and  Short 

Separation  Bubbles  on  the  Pressure 
Distribution  Over  an  Airfoil  at 
High  Reynolds  Number  (Reference 
41). 
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Figure  9.  The  Effect  of  Short  and  Long  Bubbles  on  Airfoil  Performance  at 
High  Reynolds  Number  (Reference  42). 
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Figure  10.  Airfoil  Section  Lift  and  Drag 
Coefficients  versus  Angle  of 
Attack  Showing  the  "High  C£^ax"  or 
"Clockwise"  Hysteresis. 


Figure  11.  Airfoil  Section  Lift  and  Drag 
Coefficients  versus  Angle  of 
Attack  Showing  the  "Low  C^x"  or 
"Counterclockwi se"  Hysteresis . 


Figure  12.  Sketch  of  Laminar  Separation  and 
Turbulent  Reattachment  Near  the 
Leading  Edge  of  an  Airfoil. 
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Figure  13.  Typical  Surface  Pressure 

Distribution  With  and  Without  a 
Separation  Bubble. 


Figure  14.  Successive  Locations  of  Wave  Fronts  in  the  Later  Stages  of 
Transition  (Taken  from  Reference  89). 
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b)  Plan  View 


Figure  15.  University  of  Notre  Dame 

Photographs  of  Smoke  Filaments 
Generated  by  the  Smoke-wire 
Technique  for  a  Cylindrical 
Leading  Edge-Constant  Thickness 
Model  at  6°  Angle  of  Attack  and  Rc 
<  100,000  (Air  flow  from  right  to 
left) 
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Figure  17.  Airfoil  Section  Design  Procedure  (Reference  4). 
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Figure  18.  Comparision  of  Theoretical  and  Experimental  Section 

Characteristics  of  the  EPPLER  387  Airfoil  (Reference  120). 
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Figure  18.  Continued 
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Figure  18.  Continued 
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Figure  19.  Comparison  of  Theoretical  and  Experimental  Section  Characteristics 
of  NASA  NLF ( 1 ) -04 1 6  Airfoil  for  Reynolds  Number  Equal  to  2  x  10® 
(Reference  121). 


Figure  23.  Plot  of  Computed  Streamlines  for  the  Miley  Airfoil  at  a  =  0°, 
Rc  =  100,000,  and  M  =  0.5  (Reference  123). 


Figure  24.  Pressure  Coefficient  Versus  Position  Along  the  Chord  of  the  Miley 
Airfoil  (Reference  123). 
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Figure  25.  Section  Lift  and  Drag  Coefficients  Versus  Angle  of  Attack: 

Corrected  for  Blockage  Rc  *  100,000;  One  Flow  Restrictor 
(Reference  140). 


Figure  26.  Plot  of  Computed  Streamlines  for  the  Wortinann  FX  63-137  Airfoil  at 
a  =  0°,  Rc  =  100,0(30,  and  M  =  0.5  (Reference  123). 


Figure  27.  Sketch  Indicating  the  Three-Dimensional  Corner  Type  of  Flow  Near 
the  Airfoil  and  the  Side  Plate. 


Figure  28. 


Unsteady  Wake 
of  Attack  and 
visual i zati on 


Produced  Behind  a  NACA  663-OIR  Airfoil  at  6°  Angle 
Rc  =  40,000  (flow  pattern  sketched  from  smoke 
photograph) . 


L  leading  edge 


Figure  29.  Smoke-Wire  Visualization  for  the  NACA  663-OI ft  Airfoil  Showing  the 
Spanwise  Flow  Structure  at  Rc  =  55,000  and  a  =  8°  (Reference  92). 


Figure  30.  Section  Profile-Drag  Coefficient  in  the  Spanwise  Direction  Using 
an  Integrating  Rake  (adapted  from  Althaus  Reference  129). 


Figure  31.  Low  Turbulence  Subsonic  Wind  Tunnel. 
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a)  Side  View  b)  Front  View 


Figure  32.  Test  Section  with  One  Piece  Airfoil  Model. 


Figure  33.  Side  View  of  Test  Section  with  Three  Piece  Airfoil  Model 


a)  NACA  663-OI8 


b)  EPPLER  61 


Figure  34.  Airfoil  Geometries 
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Figure  35.  Comparison  of  Lift  and  Drag  Data  Coefficients  for  the  fine  Piece 
and  the  Three  Piece  Airfoil  Models  at  Rc  =  61,400  (Reference 
127). 
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36.  Comparison  of  Lift  and  Drag  Data  Coefficients  for  the  One  Piece 
and  the  Three  Piece  Airfoil  Models  at  Rc  =  306,000(Reference 
127). 
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Figure  33.  Experimental  Static  Pressure 

Distribution  for  the  NACA  663-O 1 R 
Airfoil  at  R--  =  90,286  and  a  = 

7°. 


Figure  39.  Comparison  of  Lift  Coefficient 

Data  at  Rc  =  80,000  for  the  Eppler 
61  Airfoil  (Reference  60). 


Figure  40.  Lissaman  7769  and  Mi  1 ey  M06-13-12S  Airfoil  Geometries. 


Figure  41.  Test  Section  Showing  the  Turbulence  Screen  and  Flow  Restrictor 
Locations  (airfoil  side  plates  not  shown). 
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Figure  42.  Turbulence  Intensity  versus  Tunnel  Velocity  With  and  Without  Flow 
Restrictors  and  Screen. 
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Figure  43.  Sound  Pressure  Level  versus  Tunnel  Velocity  With  and  Without  FI 
Restrictors. 
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Figure  44.  Power  Spectral  Density  versus  Frequency  at  x  =  68  cm,  II  =  9  m/s 
and  No  Flow  Restrictor  From  Sound  Measurements. 


Figure  45.  Power  Spectral  Density  versus  Frequency  at  x  =  68  cm,  II  =  9  m/s, 
and  No  Flow  Restrictor  From  Hot-Wire  Anemometer. 


Figure  46.  Power  Spectral  Density  versus  Frequency  at  x  =  68  cm,  U  =  9  m/s 
and  One  Flow  Restrictor  From  Sound  Measurements. 


Figure  47.  Power  Spectral  Density  versus  Frequency  at  x  =  68  cm,  U  =  9  m/s 
and  One  Flow  Restrictor  From  Hot-Wire  Anemometer. 


Figure  48.  Power  Spectral  Density  versus  Frequency  at  x  =  68  cm,  II  =  9  m/s.  No 
Flow  Restrictor  and  One  7 .09  meshes/cm  Screen  From  Hot-Wire 
Anemometer. 


Figure  49.  Power  Spectral  Density  versus  Frequency  at  x  =  68  cm,  U  -  9  m/s, 
One  Flow  Restrictor  and  One  7.09  meshes/cm  Screen  From  Hot-Wire 
Anemometer. 
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Figure  50.  Lift  and  Drag  Coefficients  versus  Angle  of  Attack  of  the  Smooth 
Lissaman  Airfoil  with  No  Screen  or  Flow  Restrictor  (Hysteresis) 


h)  Decreasing  Angle  of  Attack 


Figure'52.  Lift  and  Drag  Coefficients  versus  Angle  of  Attack  for  the  Smooth 
Lissaman  Airfoil  with  No  Screen  and  One  Flow  Restrictor 
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Figure  53.  Lift  and  Drag  Coefficients  versus 
Angle  of  Attack  for  the  Smooth 
Lissaman  Airfoil  with  One  7.09 
meshes/cm  Screen  and  No  Flow 
Restrictor 
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Figure  54.  Lift  and  Drag  Coefficients  versus 
Ang’e  of  Attack  for  the  Lissaman 
Airfoil  with  Tape  Trip  at  1.1% 
Chord  and  No  Screen  or  Flow 
Restrictor 
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Figure  55.  Section  Lift  and  Profile  Drag 
Coefficients  versus  Angle  of 
Attack  at  Rc  =  150,000  for  the 
Smooth  Mi  ley  Airfoil  (Standard 
Tunnel  Configuration). 
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Figure  56.  Section  Lift  and  Drag  Coefficients 
versus  Angle  of  Attack  for  Rc  = 
200,000  for  the  Smooth  Mi  ley 
Airfoil  (Standard  Tunnel 
Configuration) . 
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a)  a  =  0°  increasing 


a  =  15°  increas^g 


c)  a  =  13°  increasing 


f)  a  =  13°  decreasing 


Figure  57.  Smoke  Visualization  Photographs  of  the  Smooth  Miley  Airfoil  at  Rc  = 
150,000  (Standard  Tunnel  Configuration). 


Figure  58.  Pressure  Coefficient  versus  Position  Along  the  Chord  of  the 
Two-Dimensional  Wortmann  FX  63-137  Airfoil  Models  Showing  the 
Influence  of  Pressure  Tap  Orientation  at  Rc  =  80,000  and  a  =  8°. 


Figure  59.  Pressure  Coefficient  versus  Position  Along  the  Chord  of  the 

Two-Dimensional  Wortmann  FX  63-137  Airfoil  Models  Showing  the 
Influence  of  Pressure  Tap  Orientation  at  Rc  =  200,000  and  a  =  16°. 
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Figure  60.  Pressure  Coefficient  versus  Position  Along  the  Chord  of  the 

Two-Dimensional  Wortmann  FX  63-137  Smooth  Airfoil  Model  at  Rc  = 
100,000  and  „  =  16°. 
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Figure  61.  Pressure  Coefficient  versus  Position  Along  the  Chord  of  the 

Two-Dimensional  Wortmann  FX  63-137  Airfoil  Showing  the  Influence 
of  Grit  Roughness  at  Rc  -  100,000  and  a  =  16°. 
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Figure  62.  The  Quarter-Chord  Moment  Coefficient  versus  Angle  of  Attack  and 
the  Lift/Orag  Polar  for  the  Smooth  Two-Dimensional  Wortmann  FX 
63-137  Airfoil  at  Rc  =  100,000. 
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Figure  63.  The  Quarter-Chord  Momentum  Coefficient  versus  Angle  of  Attack  and 
Lift/Orag  Polar  for  Two-Dimensional  Wortmann  FX  63-137  Airfoil 
with  Sawtooth  Roughness  at  Rc  =  100,000. 


a)  Section  Lift  and  Profile  drag  Coefficients  versus  Angle  of 
Attack 
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h)  Quarter-Chord  Moment  Coefficient  versus  Angle  of  Attack 
and  the  lift/Orag  Polar 


Figure  64.  Two-Dimensional  Section  Lift,  Drag,  Quarter-Chord  Moment 

Coefficients  versus  Angle  of  Attack  and  Lift/Orag  Polar  for  the 
Wortmann  FX  63-137  Smooth  Airfoil  at  Rc  =  80,000. 
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a)  Section  Lift  and  Profile  Drag  Coefficients  versus  Angle  of 
Attack 


b)  Quarter-Chord  Moment  Coefficient  versus  Angle  of  Attack 
and  the  Lift/Drag  Polar 


Figure  65.  Two-Dimensional  Section  Lift,  Drag,  Quarter-Chord  Moment 

Coefficients  versus  Angle  of  Attack  and  Lift/Drag  Polar  for  the 
Wortmann  FX  63-137  Smooth  Airfoil  at  Rc  *  100,000. 
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a)  Section  Lift  and  Profile  Drag  Coefficients  versus  Angle  of 
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b)  Quarter-Chord  Moment  Coefficient  versus  Angle  of  Attack 
and  the  Lift/Orag  Polar 


Figure  66.  Two-Dimensional  Section  Lift,  Drag,  Quarter-C.iord  Moment 

Coefficients  versus  Angle  of  Attack  and  Lift/Drag  Polar  for  the 
Wortmann  FX  63-137  Smooth  Airfoil  at  Rc  =  200,000. 
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a)  Lift  and  Profile  Drag  Coefficients  versus  Angle  of  Attack 
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b)  Quarter-Chord  Moment  Coefficient  versus  Angle  of  Attack 
and  the  Lift/Drag  Polar 


Figure  67.  Rectangular  Finite  Wing  Lift,  Drag,  Quarter-Chord  Moment 

Coefficients  versus  Angle  of  Attack  and  Lift/Drag  Polar  for 
Semi-Span  Aspect  Ratio  of  2.7  and  Rc  =  100,000. 
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Figure  68.  The  Effect  of  Semi-Span  Aspect  Ratio  on  the  Rectangular  Wing  Lift 
and  Drag  Coefficients  versus  Angle  of  Attack  for  Rr  =  200,000. 
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igure  69.  The  Effect  of  Chord  Reynolds  Number  on  Lift  and  Drag  Coefficients 
versus  Angle  of  Attack  for  a  Rectangular  Wing  with  Semi-Span 
Aspect  Ratio  of  2.0. 
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Figure  71.  Minimum  Drag  Coefficient  versus  Chord  Reynolds  Number  for  the 
FX  63-137  Airfoil  and  Rectangular  Wing. 
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Figure  72.  Pressure  Di  stributions  at  Seven  Positions  (y)  Along  Semi-Span  (b) 
and  the  Two-Dimensional  Pressure  Distributions  for  Rc  =  80,000, 

AR  =  2.0  and  a  =  18°. 
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Figure  73.  Pressure  Distributions  at  Seven  Positions  (y)  Along  Semi-Span  (b) 
and  the  Two-Dimensional  Pressure  Distributions  for  Rc  =  200,000, 
AR  =  2.0  and  a  =  18°. 
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Figure  74.  Finite  Wing  Lift  Coefficient  versus  Non-Dimensional  Distance  (y) 
Along  Semi-Span  (b}  in  Percent  for  Various  Angles  of  Attack  for 
Rc  =  200,000  and  AR  =  2.0. 
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Figure  75.  Comparison  of  Theoretical  Design  and  Experimental  Pressure 

Distribution  for  the  Miley  Airfoil  at  a  =  11°  and  Rc  =  600,000. 


Figure  76.  An  Example  of  a  Solar  High  Altitude  Powered  Platform  in  the 

Daytime  (Wing  tips  fold  down  to  increase  aspect  ratio  at  night) 
Designed  by  the  Lockheed  Missiles  and  Space  Company 
(Reference  162). 


Figure  77.  Mini-Sniffer  RPV  (Reference  5). 


a)  Horizontal  Axis  Propeller-Type  b)  Verticle  Axis  Darrieus-Type 

Figure  78.  Wind  Turbine  Configurations 


Figure  79.  Ratio  of  Blade  Tip  Speed  to  Wind  Speed  for  Various  Wind  Turbines. 
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Figure  80.  Horizonal  Axis  Mind  Turbine  (HAWT) 
Blade  Element  Velocities  and 
Forces. 
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Figure  81.  Verticle  Axis  Wind  Turbine  (VAWT) 
Blade  Element  Velocities  and 
Force. 
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Figure  82.  Thrust  Coefficient  versus  Angle  of  Attack  for  the  Wortmann 

FX  63-137  Airfoil  at  Rc  =  125,000  and  the  NACA  0015  Airfoil  at 
Rc  =  350,000. 
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